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ABSTRACT 

- T— I I 

^ ■ A description of the Fortran program HECTOR for a variety of semi-analytical calculations of radiative 
. QED, QCD, and electroweak corrections to the double-differential cross sections of MC and CC 
deep inelastic charged lepton proton (or lepton deuteron) scattering is presented. HECTOR originates 
from the substantially improved and extended earlier programs HELIOS and TERAD91. It is mainly 
intended for applications at HERA or LEP(8)LHC, but may be used also for ^.N scattering in fixed 
target experiments. The QED corrections may be calculated in different sets of variables: leptonic, 
hadronic, mixed, Jaquet-Blondel, double angle etc. Besides the leading logarithmic approximation 
up to order 0{a'^), exact 0{a) corrections and inclusive soft photon exponentiation are taken into 
account. The photoproduction region is also covered. 
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Hector, son of the Troyan King Priam. 
Troyan war leader. Homer describes Hec- 
tor to be open, frank, brave, cheerful in ad- 
versity, and tenderly compassionate in the 
Ilyad Drawing by J. Flaxman 0. 

PROGRAM SUMMARY 

Title of program: HECTOR 

Version: 1.00 November 1995 

Catalogue number: 

Program obtainable from: 

[http:/ /www. ith.de/ theory] or on request from 

e-maih hector@ifh.de 

Licensing provisions: non 

Computers: all 

Operating system: all 

Program language: FORTRAN-77 

Memory required to execute with typical data: 

4.2 Mb. 

No. of bits per word: 64; 

the branch TERADLOW uses 128 bits per word 

for numerical quadric precision 

No. of lines in distributed program: 20.000 



Other programs called: 
libraries KERNLIB [1]; 
FFREAD - part of PACKLIB [2]; 
DIZET - an electroweak library [3]; 
optionally possible: PDFLIB [4]. 
External files needed: 

HECTOR . IN? - input cards to be read by the 
FFREAD package. 

Keywords: deep inelastic scattering, QED, 
electroweak and QCD radiative corrections, 
structure functions, higher order corrections. 
Nature of physical problem: 
First and higher order QED radiative cor- 
rections to the lepton nucleon deep inelastic 
scattering; virtual electroweak and QCD cor- 
rections to the process; QCD corrections to 
structure functions. Method of solution: 
Numerical integration of analytical formulae. 
Restrictions on complexity of the problem: 
Only selective experimental cuts are possible. 
Results for full calculations in order 0{a) are 
not available for all possible kinematical vari- 
ables. 

Typical running time: 

The running time strongly depends upon the 
options used. One finds e.g.: LLA, lep- 
tonic variables, AfC, 25 points: about 25 
sec. (CPU time), full 0{a) + 0{a^) LLA 
-|- soft exponentiation + electroweak non- 
running couplings, leptonic variables, AfC, 25 
points: about 600 sec. 
References: 

[1] CERN Program Library Z 001. 

[2] R. Brun, et al., FFREAD User 
Guide and Reference Manual, CERN 
DD/US/71, CERN Program Library I 
302, February 1987. 
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LONG WRITE-UP 



1 Introduction 

HECTOR is a pro. ram to calculate radiative corrections (RC) to the charged lepton 
nucleon neutral and charged current deep inelastic scattering (DIS) l^N I'X . This requires 
to take into account the full set of virtual electroweak (EW) corrections, QCD corrections, 
and the complete QED corrections in lowest order, as well as the leading higher order QED 
corrections. The QCD radiative corrections are implemented up to next-to-leading order 
(NLO). They depend on the factorization scheme used. The 0{a) QED radiative corrections 
due to radiation from the lepton lines are known to be large and can reach the order of 100% 
in some regions of phase space. This makes, at least in parts of the phase space, the inclusion 
of higher order corrections indispensable. The aim of the project is to obtain a description at 
the percent level. 

The program originates from two different previous codes: HELIOS and TERAD91 @]. 
HELIOS calculates QED RC in the leading logarithmic approximation (LLA) in the 0{a) order. 
The program is based on refs. ||5[-0|. The 0{a^L'^) corrections and soft photon exponentiation 
have been included recently [H. TERAD91 calculates complete QED and the EW corrections 



in the 0{a) order. It is based on a series of papers 0-|14]. The current Fortran program 
HECTOR combines the two methods. 



2 Basic notations 
2.1 Born cross sections 

HECTOR allows the calculation of neutral current (A/'C) and charged current (CC) Born cross 
sections. 

The MC cross section may be calculated in several ways and approximations: in the ultra- 
relativistic approximation or exact in all masses; neglecting the longitudinal structure functions 
or not; without or with electroweak one loop and higher order corrections; the latter case 
is called improved Born cross section. The Born cross section is calculated in subroutine 
sigbrn.f : 



NC 



dxdy 



27ra^(Q^)g 



2(l-.)-2x.-+ 1-2 



+ xy{2-y)J^^{x,Q' 



Ax' 



,M2 



(2.1) 



We also introduce a slightly modified notion for later use: 



dxdy 



27ia\Q^)S 



+ 



1 



j-,Sf{y,Q') + -^,Sf{y,Q') 



2^2 



S^iy,Q')J',ix,Q' 



1 + 4x^ — 



g2 / 1 + i? 



(2.2) 
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with 



Sfiy,Q') = Q^-2m^ 
Sfiy,Q') = 2[il -y)S'-M'Q'] 
Sf{y,Q') = 2Q\2-y)S. 



(2.3) 
(2.4) 
(2.5) 



Here, M and m are the proton and lepton masses, s = —{ki + pi)"^ = S + + ^ AE^Ep, 
X and y denote the Bjorken variables, and = —Q^ is the four momentum transfer squared. 
These variables will be discussed in more detail in appendix A. The generalized structure 
functions are denoted by J-'i{x, Q^) and R is the ratio of the cross sections with virtual exchange 
of longitudinal and transverse photons, respectively: 



i?(x,Q2) 



ax 



1 + Ax- 



g2 J 2xJ-i(x,Q2 



(2.6) 



These relations are thoroughly used by HECTOR. We mention the ultra-relativistic approxima- 
tion of the Born cross section: 



dxdy 



Y^T2[x, Q2) + Y-xT^{x, g2) - y'TL^x, Q') 



with 



Y^{y) = l±{l-yf. 
For the CC reaction, a similar cross section formula holds: 



dxdy 



GlS 
An 



1 2 



F+W2(x, g2) + Y-Wsix, Q^) - y^WL{x, 



(2.7) 



(2. 



(2.9) 



The generalized structure functions for the NC case (JF) and for the CC case (W) are 
defined in section |2.2.1| , while the relations of the basic input parameters follow in the next 
section, and the description of the running of a in section p.l.2| . 



In a lowest order description one may assume the validity of Callan-Gross relations [|T5 
for 2 and Wi 2'- 



Tl{x,Q'') = ^2(a;,Q') -2x^i(x,Q2) : 
Wl{x,Q^) = W2(a;,Q2) _2xWi(x,Q2) 



0, 



(2.10) 



2.1.1 Input parameters 

HECTOR takes into account that not all the physical parameters are independent of each other 
in the Standard Model. There are three sets of input values foreseen. The default set is chosen 
in file hecset.f with flag setting IM0MS=1. For this, the input is expressed by the following 
parameters: 



fine structure constant a = 1/137.0359895; its running is discussed in section 2.1.2 
muon decay constant = 1.166388 ■ 10^^ GeV~^; 
Z boson mass Mz = 91.175 GeV; 



7 



• fermion masses, including the top quark mass rrit = 180 GeV; 

• Higgs boson mass Mh, for which the default value of 300 GeV is assumed; 

• QCD parameter A^, which is chosen in file alpqcd . f in accordance with the pre-selected 
PDF set (set of parton distribution functions). 

The default values of constants are set in file hecset.f. The strong coupling constant is 



calculated from as described in section |2.3.1 



Several additional parameters may be calculated from the above input. The most impor- 
tant ones are the W boson mass: 



Mw = Mz 
and the weak mixing angle: 



1 1 
2^2 



74.562 



1 



1 - Ar 



1/2 ^ 1/2 



(2.11) 



sin^^ 



w 



Ml 



(2.12) 



Details on the definition of Ar within HECTOR may be found in [0, [16[. It depends on the 
chosen electroweak parameters and on the fermion masses. 

If one wants to use the W mass as an input parameter instead of Mz-, the default setting 

is: 

• Mw = 80.22 GeV. 
Then, one has to determine Mz by iteratively solving the equation: 



M. 



w 



/74.562 
V M, 



w 



.-1/2 



1 - Ar 



(2.13) 



This is done when using IM0MS=2. 

A third possibility, IM0MS=3, is to take both Mz and M^r together with a as an input and 



to calculate G^: 



G,, 



\/2M2 sin^ Owil 



Ar) 



(2.14) 



For any setting of the fiag IMOMS, it is also possible to use an effective weak mixing angleQ 
instead of the definition ( |2.12| ). This is done by choosing the fiag IWEA=0, which switches off 
the calculation of virtual weak corrections. The weak mixing angle as defined in (|2.12|) is not 
used for the calculation of weak neutral couplings in this case: 



sin^ e 



w 



In file hecset .f a default value is given: 



sin^r^^ 



w 



sin^^^. 



0.232. 



(2.15) 



(2.16) 



1 For more details, see the discussion in on the definition of sin'^ 6*^. 
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2.1.2 The running electromagnetic coupling a(Q^) 



The running of the QED coupling is taken into account if IVPL=1. Then, instead of a, the 
value 



is used. The correction consists of the pieces: 



a 



1- Aa 



Aa 



Aai + Aaudcsb + Aat 



(2.17) 



f2.18) 



corresponding to the charged leptons, the lighter quarks (u-b), and the top quark. 

The value of Aa is calculated in function XFOTFl of the package DIZET. The contribution 
due to leptons is: 



Aai 



Y: Q)NfAFf{Q% 

Am) 1 / 2m2\ (3f + l\ 




1 + 



Am) 



(2.19) 



(2.20) 



f2.211 



The color factor, Nf, and the charge squared, Q), are unity for leptons. In the weak library, 
AFf is calculated by the function XI3: AFj = 2 XI3. The contribution from the light quarks 
has been parametrised in two different ways. It may be calculated by ( |2.19| ) setting IHVP=2. 
In this case, we use effective quark masses |18|: m„ = = 0.041, = 0.15, mc = 
1.5, = 4.5 GeV. The preferred, default approach (with IHVP=1) uses a parameterization of 
the hadronic vacuum polarization [|l^, which is contained in file hadrS.f. An older parame- 
terization may also be accessed choosing IHVP=3. 
Finally, the t quark corrections are: 



Aat = QtA/'tAFt(g2) + Aa2i°°P'""s_ 



f2.22) 



The last term contains higher order corrections and is calculated by the functions ALQCDS 
(approximate) or ALQCD (exact): 



AQ,21oop,aOs 



7^e^ 



VF 




(2.23) 



UY^{Q'^) is a two loop self energy function [21, 22|. For IQCD=0 (default value), this tiny 
correction is neglected. For IQCD=1 approximate and for IQCD=2 exact calculations according 
to ( |2.23| ) are made. The first one is a crude approximation only, the second, on the other 
hand, results in a very time consuming computation. 



2.2 Structure functions 
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2.2.1 Generalized structure functions 



The generalized structure functions describe the electroweak interactions of leptons with beam 
particle charge Qi and lepton beam polarization ^ via the exchange of a photon, Z boson, or 
W boson with unpolarized nucleons. 

The generalized structure functions are calculated for the MC cross sections in file genstf . f 
or, for the low branch TERADLOW, in file gsflow.f, and for the CC cross sections in file 
gccstf .f. The exact treatment of virtual electroweak corrections (for the choice IWEAK=1) 
will be described in section |2.2.3| . Here we assume IWEAK=0. Then, the generalized MC 
structure functions J^i{x, Q"^) are expressed in terms of the usual hadronic structure functions 
Fj, Gj, Hi by the following equations: 



^1,2(X,Q') 

xJ^3(x,Q^) 



+ 4 + al + 2Xveae) x\Q^)Hi,2{x, Q'), 
-2 sign(go|iai {a, + \Ve)xiQ^)xGsix,Q^) 



(2.24) 



+ 



(2.25) 



with Qe = —1 and 



A 



Op 



1 
1 



4sin2 



(2.26) 
(2.27) 
(2.28) 



and 



G,, 



The structure functions Fi^2, Gi^2, -f^i,2 and Gs, ifs are defined in the next section. 
The charged current generalized structure functions for IWEAK=0 are: 



W2(X,Q2) 

W3(x,g2) 



1 + A 



-sign{Qi)^-^xW^'{x,Q^). 



(2.29) 



(2.30) 
(2.31) 



They well defined in the subsequent section. 

For the case of an exact inclusion of the electroweak virtual corrections, we refer to sec- 



tion |2.2.3| . There we will see that the notion of the generalized structure functions holds 
further on, but that of hadronic structure functions will be lost. 



The photon exchange parts of the generalized structure functions have to vanish at low . 
This is ensured by non-perturbative modifications of the corresponding structure functions as 
described in section 
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2.2.2 Structure functions 

There are several ways to choose structure functions in HECTOR. This will be described in detail 
in section ||. 

One option is to derive structure functions from parton distribution functions (PDF). The 
different structure functions are in a lowest order approach^ 



F,{x,Q') = xJ2\Q,\'[q{x,Q')+q{x,Q% (2.32) 

G2ix,Q') = xJ2\Q,\v,[qix,Q') + qix,Q% (2.33) 

H,{x,Q') = xJ2\{vl + al)[q{x,Q')+q{x,Q% (2.34) 

xGsix,Q^) = xJ2\Q,\a,[qix,Q^)-qix,Q% (2.35) 
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xHs{x,Q') = xJ2lv,a,[q{x,Q')-q{x,Q% (2.36) 

9 ^ 

where q{x,Q^) {q{x,Q^)) denote the quark (antiquark) distributions. The neutral current 
couplings are in our notations: 

Vf = l-4|Qy|sin2^^, (2.37) 

ttf = 1. (2.38) 

Here, 9^ is the effective weak mixing angle introduced in section and SelQel = Qe = — 1 
is the electron charge. 

The CC structure functions are: 

W^ix,Q^) = 2xY.[diix,Q^)+Uiix,Q^)], (2.39) 

i 

W2-ix,Q^) = 2xY.[uiix,Q^)+diix,Q^)], (2.40) 

i 

xW+{x,Q^) = 2xJ2[d,{x,Q')-u,{x,Q^)], (2.41) 

i 

xW3'ix,Q^) = 2xJ2[uiix,Q^)-di{x,Q^)], (2.42) 

i 

with Ui and di denoting the up {u, c, t) and down quark (d, s, h) densities, respectively. 

The structure functions Si{x^Q^) = Fi{x,Q'^),Gi{x,Q'^),Hi{x,Q'^),Wi{x,Q'^) are calcu- 
lated in HECTOR by 

2xSi = S2-Sl, (2.43) 

and analogously for the other distributions. The longitudinal structure functions vanish in 
lowest order of QCD. 

^NLO corrections are delt with in section b.3l 
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2.2.3 Generalized structure functions for improved Born cross sections 

For the introduction of the generahzed structure functions in the presence of virtual weak 
interactions, we use a slightly different notation compared to that of section |2.2.1| in order to 
make the relation between HECTOR and and |Ty] sufficiently transparent. To these references 
(and for later improvements also to we refer for all the details which cannot be described 
here. 

In the MC case the structure functions are given by: 

J-f-^(x,g2) = il JCiB,p)ViB,p)q,ix,Q'), (2.44) 

B='r,I,Z i=l q=Q,Q 

Tl^{x,Q') = pJCiB,p)AiB,p)q,ix,Q'). (2.45) 

B=^,I,Z i=l q=Q,Q 

Here, index i denotes the quark generation. B = ^ stands for the pure photon exchange 
contribution, Z for the Z exchange and / for their interference. The parameter p depends on 
the relative charge assignments of the beam particle and the struck quark: 



p = I for : e Qi, e^Qi, 

p = -1 for : e+Qi, e~Qi. 



(2.46) 



The function /C contains overall normalizations: 



/C(7,P) = (QeQg?, (2.47) 
/C(/,p) = 2\Q,QQ\xmPz{p), (2.48) 
)C{Z,p) = IxmPzip)]'- (2.49) 

The function xiQ"^) has been defined in ( p.29| ). The real-valued weak form factor pz{p) was 
introduced as p^-j (p) in . The form factors depend on a variety of variables: 

pzi+l) = pz{xS,Q'^,xU]e,q), 

Pz(-l) = PzixU,Q^,xS;e,q), (2.50) 
U = Q'^-S. (2.51) 

An explicit expression for pz{S, Q^, U]a,b) in 0{a) including the ZZ, WW box terms, which 
are crucial here, is given by eq. (A.l) in 0. 

The dependence on the weak mixing angle and its radiative corrections, but also that on 
a longitudinal beam polarization ^, is contained in the functions V and A: 



V(7,p) 


= 1, 


A{i,p) 


= 0, 


V{I,p) 




A{I,p) 


= {ae + Xve) a, 


V{Z,p) 


( 2 2\ 2 
= \Ve + «e j S 


AiZ,p) 


= 2 OeOg {VeVq 
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'Vg + vl + 2Aae [Veal + VegVg j , 



(2.52) 
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The vector and axial vector couplings contain finite radiative corrections from the real- valued 
weak form factors Haip)'- 

de = 1, 

Qg = 1, 

Ve = 1 - 4: Sin^ 9w\Qe\f^e{p), 

Vq = 1 - isin^ 9w\Qg\Kg{p), 

Veq = Ve+Vg-l + 16sin^9w\QeQq\l^eq{p)- {2.53) 

The weak form factors Ka{p), a = e,q, eg, are defined in eqs. (A. 2) and (A. 3) of 0: 

Ka{+1) = Ka{xS,Q'^,xU;e,q), 

/ta(-l) = KaixU,Q^,xS;e,q). (2.54) 

Some universal higher order corrections to pz and Ha may be found in |T6|. The form factors 
are calculated in DIZET, file rokap . f , and transferred to other parts of the package as elements 
of a vector XFF. 

An approximation to the exact electroweak expressions in terms of an effective weak mixing 
angle may be introduced as follows: 



P 



"eq 



Kg ~ Kg 



1, 



together with the definition 



sin^ 9^ = K sin^ 9w- 



w- 



(2.55) 



(2.56) 



(cf. sections [2.2.1] , and \l.2.2\) . In HECTOR a numerical value may be chosen for sin^^^y. 

This value may be identified by that used in the description of the Z resonance physics at 
LEP 1: 



sin^ 9^ = Kf. sin^ 9wi 



(2.57) 



where the form factor /te describes the decay of the Z boson into two charged leptons. Thus, 
a certain scale has been chosen (Q^ = — M|), there are no contributions from box diagrams 
included etc. For details see [|T6|, ^ p^ . 

For the e~p CC scattering, the corresponding formulae are: 



xWf^{x,Q^) 



whereas, for the e^p scattering: 



xWt^^\ 



-X 



i=l 
3 



-X 



[pU^)u,{x,Q') - pU-m{x,Q' 



i=l 



-X 



i=l 
3 



-X 



J2[pU-^)di{x,Q') - pla)uiix,Q' 



1=1 



(2.58) 
(2.59) 

(2.60) 
(2.61) 
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For IM0MS=1,2, the real- valued weak form factors pc are defined as follows: 



pc'(+l) = I + 6pw{xS,Q'^,xU]e,q), 



(2.62) 



where the arguments e, q indicate the isospin and charge contents of the scattering particles. 
For IM0MS=3, the definition is instead: 



Pc{+1) 
Pc(-l) 



[1 + 5pw{xS, Q^, xU ; e, q)] 
1 - Ar ' 

[1 + 6pw{xU, Q^, xS;e, q)] 
1 - Ar ■ 



(2.63) 



The weak correction Spw is a part of the electroweak virtual CC corrections, which in- 
clude also electromagnetic contributions. The latter have to be combined with real photonic 
bremsstrahlung in order to cancel infrared divergent intermediary terms. As is well-known, 
the separation of a weak and an electromagnetic part of the virtual corrections. 



p^y = 1 6pw + Sp 



QED 



(2.64) 



is not unique for the charged current. For practical purposes, a recipe was proposed in [|lO[| , 
which allows a reasonable definition of both parts. In this reference, further details and explicit 
expressions for the form factors may be found in appendix A. 

For all the weak form factors p and k, the masses of the light leptons and quarks may be 
neglected while the t quark mass has to be taken into account. The exception is Aa, which 



depends on the light fermion masses logarit 
which is discussed later, terms of the order In 
we again refer to the quoted literature. 



imically. To the QED part of the corrections. 



S, U}/m] 



will also contribute. For details 



2.3 Parton distributions and QCD corrections to structure func- 
tions 

Nucleons and hadrons will be assumed to be built out of quarks. The quark distributions are 
labelled as follows: 



{u, d, s, c, b, t), 
(m, d, s, c,b,t), 



for i = 1, . . . , 6, 
for i = 1 6. 



(2.65) 
(2.66) 



These quark distributions contain both the valence and sea quark contributions. The explicit 
expression for xW-^ , which was introduced in (|2.42| ) may serve as an example: 



xW^{x,Q^' 



Here, we use the following notations for the parton densities: 



X 



qi{x,Q'^) ± qj{x,Q^ 



(2.67) 



(2.68) 



The combinations qfji^x^Q"^) are defined in eq. ( p.80| ). The quark densities qi are functions of 
X and Q^. In HECTOR, there are three approaches to the parton distributions available: 
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• Leading order (LO) distributions; 

• Next to leading order (NLO) distributions in the MS scheme; 

• Next to leading order distributions in the DIS scheme. 

The NLO corrections are only active if the flag IWEA = is set. Otherwise only LO distribu- 
tions are accessed. 

2.3.1 The strong coupling constant 

For the renormalization of the strong coupling constant, a^, the MS scheme is choosen. The 
expressions for as in leading and next to leading order are 



An 



/5oln(QVA2 



A lnln(QVA^ 
M[ln(gVA2)]^J 



(2.69) 
(2.70) 



with 



/So = ll-giV/, 

A = 102 -yiV^. 



(2.71) 
(2.72) 



Here, A = denotes the QCD parameter. A|^ depends on the number of active flavours, 
Nf. On passing the respective flavour thresholds the value of A is changed leaving as contin- 
uous. In the case of distributions contained in PDFACT this treatment of as is provided. For 
other cases the corresponding function ALPQCD(ISSET,Q) has to be added by the user. The 



mass thresholds in the case of the distributions 25 were choosen as in 



2.3.2 Parton distributions in the MS factorization scheme 

Several types of parton densities emerge in the structure functions and were defined in ( |2.68| ): 
qtj{x,Q^), qij{x,Q^), and q{j{x,Q^). 
The distribution g,^ reads 



where 



271 



X y 



-C i- 

y '\y, 



itiiy ^ Q^) MS -itji^^Q 



'MS 



-(itjix^Q^)Ms dyCg{y,Q^ 



xG 

Cf 
1 



(2.73) 

(2.74) 
(2.75) 

(2.76) 
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with Cf = {N^ - 1)/{2N^) = 4/3 
Correspondingly q~j is: 



qr^{x,Q')=x[q,{x,Q')-q,{x,Q')]. (2.77) 



it] 

For this combination one has: 

asiQ"^) j dyx 



y y 



-q7,{x,Q')Ms£dyCM^, (2.78) 

with: 

C,iz) = Cgiz)-il + z)Cp. (2.79) 
Finally, for the longitudinal structure function, the correction reads: 

qU^,Q') = ^^J^j[^)'[2CFqt,{y,Q')m + ^[^^ 

In a series of parameterizations all active quark flavours are delt with as massless. In these 
cases the above relations (and those given in section |2.3.3|) yield the complete description. 



On the other hand, some parameterizations account for mass effects explicitely (e.g. p6| 



leading to coefficient functions different from (|2.73| , p.78| , p. 801 ) for heavy flavour contributions 



q = c,b,t. In general these contributions are process dependent. HECTOR foresees a simple 
approach for parameterizations contained in PDFACT in the present version: the light flavours 
(m, d, s) are parametrized according to ||2B[, while the heavy flavours (c, b) are taken from [^, 
which is some suitable approximation (cf. ||26|| ). 

In all other cases (including the use of PDFLIB) the user has to provide the correct mappings 
her /himself. 

The distributions [qi{x,Q^) ± qi{x,Q^)]y^ and G{y,Q'^)j^ have to be taken from one of 
the corresponding libraries. 

In the above we have chosen as the factorization scale always. 

2.3.3 Parton distributions in the DIS factorization scheme 

The + distributions ( ^.68| ) are associated with the coefficient function 

6{l-z) (2.81) 

in the DIS scheme. Thus the generalized structure functions ^^2, ^2, ^^2, VV2, are just linear 
combinations of the quark quark distributions in the DIS scheme. The generalized structure 
functions ^3,7^3, VV3, on the other hand, constitute of linear combinations of 

Qi.ix, Q') = qr^ix, Q'Us + [ ^^Cr q^.iy, Q'Us- (2.82) 

with 

C^^^{z) = -Cf{1 + z). (2.83) 

The expression for the longitudinal structure functions are identical in the MS and DIS 
schemes. 

Again the distributions [qi{x, Q'^)±qi{x, Q^)]t)1S and G{y, Q'^)bis have to be taken from one 



of the corresponding libraries. The heavy flavour distributions are delt with as in section 2.3.2 
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2.3.4 Structure functions in different schemes 



The physical expressions of the structure functions ( p.32| )-( p.36|) , ( |2.39| )- (|2.43|) in the different 
schemes are obtained by hnear combinations of the functions: 



<,MS(Dis)(^'^')' (2-84) 

%MS(Dis)(^'^')' (2-85) 

4,MS(Dis)(^'Q')' (2-86) 
replacing the corresponding leading order expressions there. 



2.4 Low modifications of structure functions and parton distri- 
butions 

In order to ensure the vanishing of the photonic parts of the generalized structure func- 
tions at low Q^, dedicated modifications of the photonic structure functions are foreseen in 
HECTOR. They are realized by a global suppression factor FVAR in the files disepm. f , genstf . f , 
prodis.f, and gsflow.f. The default setting is IVAR=0 corresponding to [^. The structure 
functions are 



J- i)n.r 



(2.87) 



= [l-exp(-aQ2)], 
Other approaches have a similar form: 



a = 3.37 GeV" 



Fyari^Q ) 



Parameterizations of the elastic form factor 



2 5 



r 



H 
4M2^ 



r 



(2.89) 

(2.90) 
(2.91) 



may be chosen (cf. e.g. [Q) as is listed in section The other structure functions are 
continued by their value at Qq to values in the range < Ql at a given value of x. 

In figure |I| we show as an example the combination of the structure function F2 from |^ 



with a low Q parameterization as proposed in |29, 30 . 
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Q^ GeV' 



Figure 1: Low interpolation for the structure function F2{x,Q'^) between F2 = 
F^^{CTEQ3) for > 5GcV^ (ISCE = 0, ISSE =1), and the low extension seleced by 
the flag IVAR=2 for < SGeV^. 
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2.5 QED corrections. Introduction 



The inclusion of photonic corrections, 

e{h) + p{pi) ^ e{k2) + X{p2) + n-f{k), (2.92) 

is crucial for a correct analysis of deep inelastic scattering. The radiative corrections may 
differ considerably comparing different sets of outer kinematical variables. If Born kinematics 
is valid, the scaling variable y, e.g., may be defined from momentum measurements of the 
leptons or from the hadrons with no difference: 

Pih 

= EliEl^. (2.94) 

Pih 



As may be seen from figure @ [|Il|], the two definitions are quite different if a photon with 
non-vanishing momentum k is emitted. The same holds true for the variables x and Q^, or y. 

In fixed target experiments mostly leptonic variables are used. At ep colliders, a lot of 
different sets of variables may be measured: the scattering angles of the lepton or the hadronic 
system, the energies of the scattered lepton or the hadrons, the transverse hadronic momentum, 
the kinetic energy of the hadrons etc. (cf. refs. Those definitions of x,y,Q'^, which 

may be used in HECTOR are introduced in appendix |^. 

The HELIOS part of HECTOR allows the calculation of leading logarithmic approximations 
(LLA) including higher order corrections to JVC and CC scattering in a large variety of different 
variable sets. 

With the TERAD part one may calculate complete 0{a) QED corrections, including soft 
photon exponentiation, in a model independent approach. The structure functions are not 




Figure 2: Configuration of the 3-momenta in the reaction e{ki)+p{pi) l'{k2)+X{p2)+'~f{k) 
in the proton rest system, pi = 0. 
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necessarily calculated in the quark parton model and the kinematics has been worked out 
exact in the proton mass and ultra-relativistic in the electron mass. 

The TERAD part of HECTOR may also be used for the calculation of complete 0{a) QED 
corrections, including soft photon exponentiation, in the quark parton model. This makes 
use of the DISEP branch. In the quark parton model, one may, in addition to the leptonic 
corrections, also calculate the lepton quark interference terms and treat the corrections from 
the quark lines. 

The complete calculations are available for both AfC and CC scattering for a series of sets 
of kinematical variables. 

For AfC scattering at very low Q^, the complete leptonic corrections described in leptonic 
variables are available from the TERADLOW branch. Here, the exact kinematics both in the 
electron and proton masses is taken into account. 

The combined use of both parts HELIOS and TERAD is also possible. Then, the complete 
0{a) corrections of the TERAD part, without its soft photon exponentiation, are used together 
with the complete higher order corrections of the HELIOS part. 

Details on the general structure of HECTOR may be seen in figure ^ and on the availability 
of variables in figure |^. 



2.6 HELIOS: The leading logarithmic approximation 

The LLA corrections are calculated by the HELIOS branch of HECTOR. They are related to 
logarithmic mass singularities, which distinguish them from the remaining corrections within 
a given order of perturbation theory. In leading logarithmic approximation, one may separate 
the terms due to initial and final state radiation. The QED corrections are composed out of 
several contributions. In HECTOR, the LLA cross section is'^: 

^2^LLA d^(T^ ^2^ini,lloop ^2^ini,2loop ^2^ini,>2,soft ^2^im,e~^e+ 
+ + + + 



dxdy dxdy dxdy dxdy dxdy dxdy 

dxdy dxdy 



The Born contribution has been discussed in section |2.1| . It is calculated in sigbrn. f and may 
or may not include weak virtual corrections depending on flag IWEAK and QCD corrections to 
the parton distributions depending on flag ISCH. For neutral current scattering both initial 
and final state LLA corrections may contribute, while for charged current scattering there is 
only initial state radiation. For the MC case, the corresponding Born diagram is shown in 
figure 0(a). The virtual corrections to order 0{a) arise from the diagram in figure ^(b) for 
the MC case. For CC scattering, it is replaced by a diagram with a photon exchange between 
incoming electron and virtual W boson. The 0{a) diagrams with real photon emission are 
shown in figure ^. 

The two first order LLA corrections have a common generic structure: 



^2 ini{fin),lloop J. f ^2^0 

— Le j dzPil^ I d{z - z^)J{x, y, Q') 



x=x,y=y,S=S 



dxdy j 



dxdy 271 J I dxdy 

^ There are also leading logarithmic corrections related to the emission of photons from the initial or 
final state quarks. They dependend on \n{Q^ /m^) iff calculated in the on-mass-shell scheme. These terms 



deserve a special treatment. One can absorb these corrections into parton distributions |[ |l^, |39) or 
fragmentation functions ||] leading to modified expressions of the scaling violations of these quantities. 
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p(l) 



(2.96) 
(2.97) 



Here, we introduced the notion 



Le = In— - 1. 

mi 



(2.98) 



The definition ( |2.98|) reproduces the soft photon terms of complete calculations in leptonic 
variables (see e.g. |Tl|). 

The expressions (|2.96| ) for initial and final state radiation differ by the scaling properties of 
the kinematical variables, the resulting Jacobian, J , and the integration bounds. In table ^ for 
the initial state corrections these parameters are listed for a variety of different sets of variables. 
Table |^ contains the same for final state corrections if they exist. For some variables, e.g. the 
Jaquet-Blondel variables, the final state is treated completely inclusive. Then, in accordance 



with the Kinoshita-Lee-Nauenberg theorem |^ there is no LLA correction. 

Furthermore, there are certain second order corrections from the initial state which may 
be non-negligible 



^2^ini,2loop 

dxdy 



a 

2^^ 



' £dzp(!''\z)\^e{z-zo)j{ 



x,y,z) 



dxdy 



x=x,y=y,S=S 



dxdy 



+ 



Q 



(Z) + Le E In 32 Pee}\^) } Ji^, ^) 



a 



f=i,q J 



dxdy 



(2.99) 



x=x,y=y,S=S 
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Figure 4: Leptonic QED hremsstrahlung diagrams 



The second order splitting functions are: 



p(2,l)( 



[21n(l-2) -ln^ + - 
l-2 [ ^ ' 2 



+ -{l + z)hiz-{l-z) 



67 



76 



{l + z)\nz + ^il-z) + ^^il-z'), 



Nc{f)Q}\p^lH 



1- z- 



2m 



The symbol ® in (p.lOO|) and ( p.l01| ) denotes the Mellin convolution: 



A{x) ®B{x) 



dxi / dx25{x — xiX2)A{xi)B{x2)- 



(2.100) 

(2.101) 
(2.102) 

(2.103) 



Qf is the fermion charge, and Nc{f) = 3 for quarks, Nc{f) = 1 for leptons, respectively. 
In (|2.102|) , we parametrize ^^^^(s, mj.) in terms of effective quark masses and use the param- 
eters m„ = 62 MeV, = 83 MeV, = 215 MeV, = 1.5 GeV, and mt, = 4.5 GeV as 
obtained in [0, [g. 

The HELIOS part of HECTOR allows to include higher order soft photon corrections into cross 
section calculations: 



^2^(>2,so/t) 



dxdy 



dzP>'\z,Q'){9{z-zo)J{x,y,z) 



dxdy 



x=x,y=y,S=S 



dxdy j 



(2.104) 



with 



D 



a 



NS 



C = 



exp 



27T 1-Z 

K (I - 



1 H Le 

271 



11 



+ 21n(l -z) 



r(i + C) 

31n[l - (a/37r)Le] . 



(2.105) 

(2.106) 
(2.107) 
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Equation ( p.l06| ) was derived in |^ already. If the charge of the final state electron is not 
recorded, also the conversion cross section: 



^2^(2,e-^e+) i 



dxdy 

with the conversion rate: 



dzP{z,Q'^;e ^ e~^)J{x,y, z] 
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dxdy 



x=x,y=y,S=S 



(2.108) 



a 
2^ 



(2.109) 



has to be included into the radiative corrections. 

The explicit expressions for the splitting functions have been taken from where addi- 
tional useful comments on the higher order corrections may be found. 

Finally, we come to the Compton correction. It has to be taken into account in the case 
of leptonic variables: 



a 



dxidyi xiS 



1 + (1-1/0' 



(2.110) 



This expression may be extracted from the complete calculation of the 0{a) corrections [§, |TT| 
and is accessed setting the flag ICMP=1. An alternative description of the Compton peak has 
been derived in (ICMP=2): 



d'a 



dxidyi 



''^D ,iz O') ^'^(67 67) 
^^^/^^^'^^^ dxdm 



(2.111) 



s=zs,x=xi/z 



with: 



d'^a{e'j — ^ 67) 
dxdyi 

D^/pixuQf) 



-5(1 -x) 



27ra^ 1 + (1 - yi) 

s l-yi 



(2.112) 



A, 7q2 _ Qi xi 



7^ ^HT'^'^ 



(2.113) 



In principle, the Compton correction has to be taken into account also in the case of mixed 
variables. But, although being logarithmically enhanced, it is not larger than the rest of the 
0{a) corrections in this case 



2.6.1 Cuts 

In the LLA approach, it is possible to apply a cut on the reconstructed electron energy by 
cutting the integration variable z. The variable ZCUT is defined by: 



[Ej{l - cos^j) + Ee'{l - cos^e')] > ZCUT. 



(2.114) 



It is also possible to reject the production of final states with a Ql and/or an invariant hadronic 
mass below some bound. See also the description of flag IHCU in section ^J, . 
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Initial state radiation 




S 




y 


zo 


J{x,y,z) 


leptonic variables 


Sz 


Qfz 


z + yi-1 

Z 


^-yi 

1 - xiyi 


yi 

z + yi-'^ 


mixed variables 


Sz 


Qfz 


yjB 

z 




1 


hadronic variables 


Sz 


Ql 


yh 

z 


yh 


1 

z 


.TB variables 


Sz 




Ujb 
z 


yp , 


1 - 2/.7B 


z-y^B 


l-iCjeil-yje) 




double angle method 


Sz 









z 


&i,y3B 


Sz 




Ujb 
z 




f - y.TB 

1 -J/jB 


S method 


Sz 


Ql 


y^ 




1 

z 


eS method 


Sz 


Qfz 


ye^z 




1 



Table 1: Scaling properties of various sets of kinematical variables for leptonic initial state 
radiation. The different parameters are defined in appendix A. 
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Final state radiation 




S 




y 




J{x,y,z) 


leptonic variables 


S 


z 


z + yi-i 

z 


1 - - xi) 


yi 

z{z + yi-l) 


mixed variables 


s 


Qf 

z 






1 

z 


S method 


s 




l-y^{l-z) 


^0 


1 

? 


eE method 


s 


z 


y T.z^ 


^0 




[i-yA^-z)f 


[1- 2/^,(1 



Table 2: Scaling properties of various sets of kinematical variables for leptonic final state 
radiation. The different parameters are defined in appendix A. 
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2.7 TERAD: Complete 0{a) leptonic corrections in the model inde- 
pendent approach for J\fC scattering 

If one wants to improve the LLA approximation of QED corrections, a natural step is to 
calculate the contributions from figures ^(b) and which lead to the leptonic corrections, 
exact in order 0{a). This is much more involved than the LLA calculation alone and has 
been done so far for a rather limited set of variables only. A comprehensive presentation of 
the calculation and of the numerical details may be found in [TTI and in references quoted 
therein. Here we restrict ourselves to a short collection of formulae, which are realized in 
HECTOR. 



Let us remind the Born cross section as given in (|2.2|) : 



2^B 



dydQ"^ 



27CQ^S_ 



(2.115) 



with Xs = S*^— 4m^M^. The factorized kinematical functions Sf are introduced in (2.3)-(2.5) 
The functions Ai are related to the generalized structure functions by: 



^2(x,Q2) 



2^i(x,Q2), 
2^^3(x,Q^ 



(2.116) 



In the TERAD part of HECTOR, the QED corrections have a structure quite similar to that 
of the Born cross section (|2.115|) : 



dS 
Si£,I) 



Jd'l-^^S{S,T), 



nXs 

^A(x,Q')5,(£:,J) 



(2.117) 
(2.118) 



Here, 8 is a set of two variables on which d'^a finally depends. In the case of leptonic, mixed, 
hadronic variables, e.g.: 

£ = {yi,QnAyh,QnAyk,Ql}- 

The X is a corresponding set of two variables to be integrated over: 

i={yh,Ql},{yi,Ql}Ayi,Qn- 

The functions Si{S,X) in ( ^.118| ) mark the difference to the Born cross section and have to 
be calculated. They are results of an one-dimensional analytical integration. In some of the 
variable sets additional integrations have been performed analytically. 

We only mention here that in the exact calculations the correct treatment of the kinemat- 
ically allowed regions of the integration variables may be crucial. A careful and exhaustive 



analysis of this may be found in appendix B of [11 



The cross sections may be represented in the following generic form: 



dS 



1 + ^ 5'yT''i£) 



+ 



d£ 



The term: 



a 

TT 



cexpon 



exp 



1 + 



a r 

71 



dS 



(2.119) 



(2.120) 
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contains the (exponentiated) virtual and soft photon part of the corrections and d'^a^^'^"^^ the 
finite part of real bremsstrahlung. Exponentiation may be taken into account with flag lEXP. 

The Born cross section is called from subroutine sigbrn.f, while the non-factorizing part 
of the QED corrections directly uses file genstf .f . The corrections may or may not include 
weak virtual corrections depending on the flag IWEAK and QCD corrections to the parton 
distributions depending on the flag ISCH. These two types of corrections are not yet thoroughly 



compatible in the present version of HECTOR as mentioned in section ^3 



2.7.1 Leptonic variables 

In leptonic variables one has: 



5 infill) - ^In^ 



1 - yi{l - xi] 



+ |ln4-Li,(i; 



1 - ym 
-2 , 



+ Li2 



i-yi 



[I - yixi)[l - yi{l - xi)] 



where 



5,nf{Si) = [in^ - l) In 



(1 - yixi)[l - yi{l - xi)] 
Here, Li2(x) denotes the Euler dilogarithmic function 



Li2(a;) 



1 , ln(l — xz) 
dz 

z 



The finite bremsstrahlung correction has a factorizing and a non-factorizing part: 



2^brems 



d'a 



dyidQi 



dyhdQl 



i=l 



Ai{xh, Ql)j7iSi{yu Qh Vh, Ql 



- A{xu Q?)^ Sr{yi, Qf) c'^^y^, Ql y,, Ql) 



The radiators Si{S,I) of the non-factorizing part are: 



-UqI + Ql) + 



Qt - 

Ql~Ql 



m\Ql-2m')^ 



1 



(l)^-(l-yO 



S2{S,X) 



^2 



+ 



M\Ql + Q',)-yn{l-yi)S' 
^ Ql{l)[{l)-yK]S' 



{Ql-Qf)VC-2 
+ (1 - yi) [(1 - m) + y,] - 2M\Ql + 2m^) 

+ 2m' [[(1) - y,] [(1 - y{) + yn] + (1)(1 - yi) 



(2.121) 



(2.122) 



(2.123) 



(2.124) 



(2.125) 
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2m 



2 B2 



2M2 

7^ 



(i)^-(i-i/,) 



(2.126) 



S 



^2 



2Ql{Ql + 2m') [(1) + (1-1/0] 



Ql-Qf 

Co J 



2(1 - m)Ql - y,{Ql + Q] 



+ (1)^ -(1-2/0 



(2.127) 



with the following abbreviations: 

B2 = {2M'Q^iQ^-Ql) + il-yi){yiQl-y,Q^)S' + S'{l)Qf{yi-y,,)} 
^ - B,{il)^-il-m)}, 

C2 = {(1 - yi)Ql - Qf [(1) - yH]fs' + Am'[im - yh)iyiQl - yhQf)S' - M^Ql - Q^Y 

^ Ci{(l)^-(l-y,)}- 

The expression for the correction factorizing off the Born term, C^^{yi, Qf ,yh,Q\), reads (exact 
in both masses m and M): 

Q2 + 2m2. 1 1 \ 2/ 



C'^'iyuQlyj^.Ql) 



Bi ^ B2 



3/2 ' ^3/2 



(2.128) 



The boundaries for the two-dimensional numerical integration in ( p.l24| ) are derived in ap- 



pendix B.2.3 of [11 



The above formulae are realized in file teradl .f . 



2.7.2 Mixed variables 

In mixed variables, the following functions have to be used in (|2.120|) : 





Xm{l - Vh) 




Xm 1 



2 m"^ 



where 



ln^-lVn[(l-l/.)(l-x„)]. 



The final formula for d'^a^^'^"^^ / dyhdQf can be written in the form: 



2^brems 



d'a 



dyndQf 



2a^S' 
Ac 



2 3 ]^ 

L„ (^Ql E [M^h, QD-T^Sjiyi, Ql yh, Ql 

i=l 



- A{x^. Ql)-^,sf{yu. QD^yu Ql Ql) 

+ /, dQlY.[A{xn.Ql)-^Sf{Qlyu,Ql 



i=l 
1 



A{x^, Qf)^,Sf{Qf, yu)C'^{Ql yn, Ql) 



(2.129) 



(2.130) 



(2.131) 
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The integration boundary for the remaining one-dimensional integration may be found 

in appendix B.3.3 of [0. The integration region must be spht into two regions / and // as fong 
as the condition Qf < y^S is fulfilled, which corresponds to < 1- If instead > 1, there 
is only the region /. The factorizing part of the cross section differs in the two kinematical 
regions: 



1 



S{Q1 - Ql) 
1 



S{Q1 - Ql) 

The following abbreviations are used: 



Qh 



l-^(Li + L2) + (L + Li + Li) 



Qh 



S(L + L,) 



(2.132) 
(2.133) 



In 



Q'l 



In 



QV 



In- 



Vh 



(2.134) 



Li = In 



Ql - VhQl 



Lo = In 



ii-yh)Ql 



I Ql Qh I \ Ql Qh \ 

The non-factorizing hard bremsstrahlung functions in the two integration regions are: 



(2.135) 



Si{QlyH,Ql) 



S',{Qly,,Ql) 



1 

S 



Q 



Q 



^Ql-Ql 



1 + ^ (L + L1-L2-1) 



Q' 



iQf 



Q'h 



Si{Qlyh,Ql) = 5. 



VkQl f, , Qt 



^Ql V ' Qt 



1 Ql 



\_QV 

: QL 



(2.136) 



Qf Ql 



Ql 



J_ Qh 
^QlV^Qf. 



i + ^j (L + L1-L2-1) 



+ 



3g| 
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Si\Qlyh,Ql) 
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+ 2(l-y.)§[ 



(2.141) 



The above formulae are realized in file teradm.f. 



2.7.3 Hadronic variables 

There exists a closed expression for the complete leptonic QED corrections in hadronic vari- 
ables to order 0{a) with soft-photon exponentiation: 



exp 



a 



dyudQl dyhdQl 
The soft photon corrections are 



TT 



5in}{yh, Ql) 



2o? 



1 



<^in/(yh,gD= (ln^-ljln(l-y,)- 



The hard bremsstrahlung corrections are taken into account to order Oipt): 



(2.142) 
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-Jln|//,Lh + 7(1 + — ) Lhi 
2 4:\ yhj 



1) 

yhJ 



-2yh^^^yh - (1 -2(l-|/h)ln7/,,ln(l-y,j) -7/ftLi2(l) 
- (2 - 3yh)U2{yh) + yn^^yh^h + y (1 - yh)Ki - y(2 - yu) ^^yh 



30 



- Li2(l) +lnyh + lnyhLh 



- In^ y,, + - yn) + 2 Inyh ln(l - yh) + 3U2{yh) 



3 7 1 

+ -^Vhl-hi + yh - - + 2 (1- 2yh) Iny^ V 



(2.143) 



with Lh = hi{Ql/m^) and Lm = U + ln[7/ft/(l - yu)]. 

Since the structure functions depend on the hadronic variables and are hardly to be in- 
tegrated over together with radiative corrections, such a compact result may be obtained in 
hadronic variables only. 

The above formulae are realized in files teradh.f and terhin.f using two different ana- 
lytical expressions. 



2.7.4 Jaquet-Blondel variables 

The net cross section is: 
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The net factorizing part is: 
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The non-factorizing part of the radiative cross section in Jaquet-Blondel variables is: 
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(2.149) 
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The integral over r is performed in HECTOR numerically since one cannot neglect the dif- 
ference between Q^^^ and Q\ - the structure functions are depending on Q\ and thus on r. 
The above formulae are realized in file teradj .f. 



2.7.5 Cuts 

In the TERAD approach it is possible to reject the production of final states with a Q\ and/or 
an invariant hadronic mass W'^ below some bound. This is realized in HECTOR for all the sets 



of kinematical variables. See also the description of flag IHCU in section 4.2 



The QED corrections in leptonic variables in the TERAD part of HECTOR allow, in addition, 
the application of cuts on the photon kinematics. Details of the formulation may be found in 
section 8.3 and appendix B.1.2 of [|ll|. In HECTOR this is done by the variable GCUT which is 
described in section 



2.8 DISEP: Complete 0{a) corrections in the quark parton model 
approach 

Complete order 0{a) corrections in the quark-parton approach are calculated in the DISEP 
part of HECTOR. They include a large variety of contributions: 

(i) Electroweak corrections, including the running QED coupling. 

(ii) Photonic corrections due to the emission of photons from the leptons together with the 
corresponding vertex corrections. 

(iii) Photonic corrections due to the emission of photons from the quark lines together with 
the corresponding vertex corrections. Their treatment deserves some care in view of 
the unknown quark masses leading to logarithmic mass singularities in the on-mass-shell 
scheme. We refer for details to the literature ||3^, ^ |3^, [13], |3^. In the branches selected 
by IDIS=1,2 these terms do not contribute. 

(iv) The interference of photon emission from leptons and quarks together with the corre- 
sponding 77, 'jZ, 'jW box corrections. 

(v) Exponentiation of soft photon corrections due to (ii). 

In the present version there are no QCD corrections taken into account in the parton distri- 
butions in this part of HECTOR. 

The additional bremsstrahlung diagrams for JVC scattering are shown in figure ^ The 
interference of them among themselves is combined with vertex corrections from the analogue 
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Figure 5: QED Bremsstrahlung off quarks 



of figure ^(b), and the interference of them with the leptonic radiation (figure ^ has to be 
combined with the 7Z and 77 boxes in the MC case and of 7!^ boxes in the CC case. 

The CC process has some specific features. The leptonic vertex correction comes from a 
diagram, where the virtual photon connects the incoming electron with the virtual W boson; 
the final state neutrino has no photonic interactions. The virtual photonic corrections from the 
quark side arise analogously. Furthermore, there are contributions from a Feynman diagram 
with real photon emission from the virtual W boson. All these contributions are taken into 
account in HECTOR. Some details of gauge invariance and other peculiarities may be found 
in 



10 



The electroweak corrections (i) are taken into account by generalized structure functions 
as is described in section |2.2.3| . 



If identical structure functions are chosen, the leptonic QED corrections (ii) agree with 
those calculated in the model independent approach of section although the underlying 
analytical expressions look quite differently. There is no access in the model independent 
approach to the corrections (iii) and (iv). 

The leading logarithmic approximation reduces the complete 0{a) corrections to terms 
proportional to ln(Q^/mg) in (ii), to ln((5^/m^) in (iii). In (iv) no mass singularities are 
contained. 

At present, in HECTOR are the following branches with QPM based complete calculations 
available: 

(i) MC and CC deep inelastic scattering in leptonic variables; 

(ii) MC deep inelastic scattering in mixed variables. 

2.8.1 He scattering in leptonic variables 

The QPM based electroweak corrections to MC scattering in leptonic variables are described 
in §. 

The virtual corrections are a part of the generalized structure functions in the improved 
Born approximation, calculated in subroutine sigbrn.f . They are described in section ^.2.3 ^ 



^Here we deviate slightly from 
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In addition, the QED corrections have to be determined. In the DISEP part, they have the 
following structure: 



^2^QED 27ra\Qf)Sxi 



dxidyi 



Qt 



Y.Y.Y.^MB,P) \V{B,p)R^{B) +pA{B,p)Rt{B)] , (2.151) 



B b 



where the sum over Q, Q extends over quarks and anti-quarks and B = {7, /, Z}. The sum over 
b covers bremsstrahlung from leptons, lepton-quark interference, and quarks {b = {e,i,q}). 

The parton distributions are accessed via the subroutine prodis.f directly. The explicit 
use of QCD corrections to the parton densities is not foreseen here. 



The parameter p and the couplings }C{B,p),V{B,p),A{B,p) are defined in section |2.2.3 
The coefficients Cb contain charges: 



Ce 
Ci 



Ql 

QeQqi 



The QED corrections are described by the functions E^''^{B): 

= Rb{B- 1, l-yi) + sl^Rb{B- y, - 1, -1), 

with 



V V 



1. 



(2.152) 
(2.153) 
(2.154) 
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The factors Rb{B; 1,1 — y) have the following structure: 

Rc{B;a,b) = -\s,iB;a,b)fg{xi,Qf)+ [\x, 

71 I Jxi 



f,{xh,Q^)Uc{B-a,b) 



+ 



rcXhfqjXh, Ql) - XlfqjXi, Qf 
Xh - Xl 



-Tc{B;a,b) 



with 



= 1/Xh, 1, 1 



for c 



e,i,q. 



(2.157) 



(2.158) 



The S, T, U are results of a two-fold analytical integrations over photonic angles. They are 
the analogues to the kinematical functions Si and shown in ( p.l25| ) to ( p.l28| ). Because 
there are many of them, we refer for explicit expressions to the original literature. The six 
independent functions are: 

Re{j;a,b), Re{r,a,b), Re{Z;a,b), Ri{-f;a,b), Ri{Z;a,b), Rq{-f;a,b). 

They are composed in the subroutines SGGiL , SGZiL , SZZiL, i=L , I , Q of HECTOR and the 
constituents S,T,U are calculated in functions FGGiL, FGZiL, FZZiL, i=L,I,Q. The others 
are expressed by the following identities: 



Ri {r,a,b) = -[Ri{'y;a,b) + Ri{Z;a,b)] 
Rq{'y;a,b) = Rq{r, a,b) = Rq{Z; a,b). 



(2.159) 
(2.160) 
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At the end of this section we should mention that the virtual corrections, by which the 
QED corrections are accompanied via the weak and electromagnetic couplings, may also be 
calculated with the use of the actual kinematical variables, i.e. folding them together with the 
functions T and U. This may be achieved by chosing IC0NV=1. The resulting computation is 
very time consuming and leads to numerical values which differ only insignificantly from those 
obtained without using this option. 

In the QPM approach in leptonic variables the virtuality of the parton densities was choosen 
by the outer variable Qf. For hard photon emission, one actually would have to perform 
integrations over weighted structure functions at reduced kinematical variables. The structure 
functions depend on the hadron kinematics, i.e. on Xh and Ql- In the DISEP approach Ql 
is to be integrated over analytically, which leads, however, to the problem that the quark 
distribution functions depend on some Q^, but not on the correct (hadronic) Ql- The same 
procedure is followed in the LLA calculations. In this sense, the DISEP approach has a reduced 
accuracy in leptonic variables compared to the model independent approach. For a high 
precision calculation it may be useful to calculate the largest corrections (ii) from TERAD and 
the remaining terms from DISEP. 



2.8.2 CC scattering in leptonic variables 



The QPM based electroweak corrections in HECTOR to CC scattering in leptonic variables are 
described in [IC]. As in the case of JVC scattering, the improved CC Born cross section is 
determined in subroutine sigbrn.f. 

The QED corrections have access to the parton distributions not through the general- 
ized structure functions of file gccstf .f but instead directly through prodis.f. The QED 
corrections are: 



dxidyi 



GlSxi 

TT 



1 + A 



l+p 



1 — p - ■ 



(2.161) 

The notations are those being introduced in sections p.2.3| and p.8.1| . The generic definition 
of the functions Rb is ( |2.157D . They are calculated in subroutines SWWL, SWWI , and SWWQ of 
HECTOR. Explicit expressions may be found in appendix B of |]10| . 



As was mentioned already in the case of AfC scattering, the QED corrections in leptonic 
variables in the DISEP approach are calculated with parton distributions having a depen- 
dence on instead of the correct variable Qf^. 



2.8.3 J\fC scattering in mixed variables 



The improved Born cross section is as described in sections 2.1 and 2.2.3 



m 



The QPM based QED corrections to ATC scattering in mixed variables have been described 
T3[1 . The structure of the corrections is the following: 



dxmdym 
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where 



+ / dz 

lym 



(^1 Um) ~\~ V^a (^i Un 



VoY+ 



+ pAoY+ 



(2.162) 



(2.163) 



Much of the notations can be found in the previous sections and section |2.2.3| . In addition, 
we use for the vector couphngs the following modifications: 



V{a,b) = K{a,b;B,p)V{B,p), 

B=^,I,Z 

K{a,b;B,p) = xbA") XB2ip) fC{B ,p), 



(2.164) 



G, Ml 
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V2 8na{Q^)aQl + Ml 



The axial couplings A(a, b) are related to the couplings A{B,p) of p.52| ) in the same manner 
as is explained above for the vector couplings. 

The functions 5*6, Sg contain factorizing soft photon corrections and the corresponding 
vertex corrections. 

The hard bremsstrahlung functions in ( p.l62|) are: 
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B a ■^mi yrri) 
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Lan+Ua 2/m) /q (^a;„) | , 

(2.165) 

with the modifications of weak neutral coupling factors; for the vector couplings: 

K = c,V{z,z), V, = CiV{l,z), Vg = c,V{l,l). (2.166) 



The bremsstrahlung axial couplings are defined analogously. For a = e,i, q, index n in ( p.l65| ) 
runs from 1 to 2, 1, 2. The functions Fan are regularized by a subtraction: 



Reg 
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m I 



(-) 



V,A 



(-) 



V,A 



Fan iz,ym)fQizX^) - Fan (1, 2/m)/Q(x 



(2.167) 



Explicit expressions for the functions S, F, L, U may be found in [1^. In HECTOR they are 
realized in subroutines SGGiM, SGZiM, SZZiM, i=L, I ,Q. 



2.9 TERADLOW: QED corrections at low 

In HECTOR, the double differential cross section of the photoproduction process is calculated in 
leptonic variables with the aid of file terlow.f. For details of the derivation of formulae and 
of the kinematics, which has to be treated exactly in both the electron and proton masses, we 



refer to 111 
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The process kinematics is described by Qf as defined from the leptonic variables and the 
invariant mass W of the system consisting both of the photon and the hadrons, 



Q^ = Qi = {h-k2 



iQi+Pif- 



(2.1Q8) 



For a fixed value of W'^ <^ S, the minimal value of Qi may become extremely small, but 
remains non-vanishing: 



[W^ - M 



212 



> 0. 



(2.169) 



This leaves the integrated deep inelastic cross section finite. The net cross section is: 



dyidQf dyidQf dyidQf dyidQf 
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The contribution from the anomalous magnetic moment is: 
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(2.172) 
(2.173) 

(2.174) 



The factorizing photonic corrections from bremsstrahlung and the usual, exact vertex correc- 
tion are: 



2^ ^ 2/3 



L^ln-:;;^ + Ll2 



- 1 
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(2.176) 
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(2.178) 
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+ /9o 



ka = kia + ^2(1 — a), 



(2.179) 
(2.180) 
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k1 



f2.1811 



Finally, the infrared finite part of the real bremsstrahlung can be written as follows: 



2^F 



a 



R 



2a^S 



dyidQf 



dMidQi y: 



i=l 



Ai{xh, Ql)j7iSi{yu Q^, Vh, Ql 
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1,2,3) are given by (^ ) - (|27[27|) , and C'^'iyi^Qlyh^Ql) is 



where Sjjyi Qf yh ,Ql) (i 
defined by ( |2.133[ ). We would like to remind that these expressions are exact in both masses 
m and M for the photon exchange contributions. 

2.9.1 Cuts 

In TERADLOW it is possible to reject the production of final states with a Q\ and/or an invariant 
hadronic mass below some bound, (cf. the description of flag IHCU in section |4.2| ). 

2.10 Other reactions 

In the LLA approach to leptonic QED corrections the Born cross section is folded with splitting 
functions at rescaled kinematic variables. This factorization may be used also for exclusive 
reactions having the same leptonic tensors as in the case of deep inelastic scattering. Examples 
for this class of reactions are e.g. heavy quark production, deep inelastic scattering including 
Z' exchange, and Higgs boson production. These processes were studied previously in this 
context g, 13, 

Sometimes one is interested even in the LLA contributions either due to initial or final 
state radiation only. For this case it is sufficient that the subsystem ('Born') cross section 
has a charged lepton line either in the initial or final state for which the LLA terms shall be 
calculated. 

In both these cases the respective corrections can be obtained setting the fiag IBRN=0. 
The user has to supply a subroutine usrbrn . f describing the subsystem process. 
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3 Structure of the Code 



GSFLOW 



HECAUX 



HECFFR 



HECOUT 



The Fortran program HECTOR consists of several subdirectories. They are shown in figure ^ 
Some of them contain several files, others only one. Some of the files contain several subroutines 
and functions. Exceptions, which are not shown in figure are: the subdirectory DIZET, 
which is used for the calculation of the electroweak radiative corrections and the subdirectory 
HECAUX with some auxiliary routines. Besides the subdirectories, the two files HECTOR. INP 
and HECTOR. OUT are added in the fiowchart. 

The subdirectories are used for the following tasks: 

CMAIN : Organization of the program fiow for combined use of the branches HELIOS and 
TERAD 

DISEP : Package for calculation of complete order 0{a) QED corrections with soft 

photon exponentiation determined in the quark parton model approach 
DIZET : Package for the calculation of weak virtual corrections 

GENSTF : Calculation of the generalized structure functions for HELIOS, TERAD, and DISEP 
Calculation of generalized structure functions for TERADLOW 
Some auxiliary functions 

Reads input data cards from file HECTOR . INP using FFREAD 
Writes output to file HECTOR . OUT, partly using FFREAD 
Setting of initialization parameters 
The main program 

Package for calculation of QED corrections in the leading logarithmic 
approximation, including higher order corrections and soft photon exponentiation 
HMAIN : Organization of the program fiow for branch HELIOS 

LOWLIB : Package for the optional modification of structure functions at low values 
PRODIS : Calculation of structure functions for electroweak corrections in all approaches 

and for QED corrections in package DISEP 
PDFLIB : Standard package of parton distribution functions 
Default package of parton distribution functions 
Calculation of the (effective) Born cross section 

Calculation of structure functions from parton distribution functions 
TERAD : Package for model independent calculation of complete leptonic order 0{a) QED 

corrections with soft photon exponentiation 
TERADLOW : Package for model independent calculation of complete leptonic order 0{a) QED 
corrections with soft photon exponentiation intended for the region of extremely 
low Q2 

TMAIN : Organization of the program flow for the branch TERAD, which may call the QED 
corrections of the packages DISEP, TERAD, or TERADLOW 



HECSET 



HECTOR 



HELIOS 



PDFACT 



SIGBRN 



STRUFC 



The user supplied subdirectories are installed containing trivial versions only: 



USRINI 



USROUT 



USRBRN 



USRSTR 



USRPDF 



File for input modiflcations, which go beyond the change of data cards foreseen 

Modifled output speciflcations 

User supplied Born cross section 

User supplied structure functions 

User supplied parton distribution functions 



The two flies in flgure || are: 
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HECTOR. INP : File with the input data card contents 
HECTOR .OUT : Standard output file 



Figure |^ is devoted to the interplay of the subdirectories GENSTF, DIZET, and LOWLIB. 
The files may be characterized as follows: 



CCOUPL : 

DIZET : 

GCCSTF 



GENSTF 



LOWLIB 



MIDSTF 



NCOUPL : 

POLYN : 

PRODIS 



STRFBS 



Calculation of charged current Born couplings or effective couplings 
Standard Model Electroweak Library 

Calculation of generalized structure functions for charged current reactions 
Calculation of generalized structure functions for neutral current reactions 
Files with different structure function modifications at low 
Interpolation between the structure function from GENSTF and the low 
modified structure function 

Calculation of neutral current Born couplings or effective couplings 
Auxiliary function for MIDSTF 

Parton distributions of the proton (used by DISEP) 
File with the Brasse/Stein modification of structure functions 
USRSTR : User supplied structure functions 

The entry labelled '3' connects the present figure with figure ^. 

We show in figure ^ how the actual structure functions are built from the parton distri- 
bution functions. The subdirectories STRUFC, PDFACT, USRPDF, PDFLIB are shown as boxes. 
Their interaction is organized by the flags ISTR, ISCH, ISSE. 

HECTOR consists of four branches as may be seen in figure P: HELIOS, TERAD, DISEP, 
TERADLOW. The first three branches of HECTOR which contain the different packages for the 
calculation of the QED corrections are shown in figure ^ The subdirectories HELIOS, TERAD, 
DISEP are symbolized by the big boxes. The settings of flags IMEA and IDSP define the choice 
of the type of corrections. HELIOS consists of one file, while inside TERAD and DISEP, the file 
structure is shown. The index i in the figure is a short-hand notation. For i=l, it is I0RD=2 
and IEPC=0, while i=2 corresponds to I0RD=2 and IEPC=1. 

The branch TERADLOW is attained by IL0W=1 together with IMEA=1 and I0PT=2. 



Figure |T0| illustrates the influence of the flags lUSR, ITVl, ITV2, IBIN, lOPT, lORD on the 
initialization of HECTOR. Shown are the subdirectories HECSET and HECTOR, together with some 
selected files (in boxes). We should comment shortly on the following of them: 



FILBIN 



HECCAL 



HECRUN 



HECTER 



USRBIN 



USRINT 



The selected lattice of kinematical points is filled 

Calls the branch of HECTOR chosen and initializes the actual calculation 
Runs the job 

The common blocks of TERAD are filled 

Optional filling of a user supplied kinematical lattice 

A user supplied integration over bins is performed 



The entries 1 and 2 are links to figure 11, where the output organization is illustrated. 
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HECOUT 
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CMAIN 
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Figure 6: Basic HECTOR flowchart. Shown are subdirectories. 
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Figure 7: Logical structure of GENSTF 
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STRUFC 



STRUFC 



ISTR 



PDFACT 
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Figure 8: Building structure functions in HECTOR 
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HELIOS 
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Figure 9: The branches of HECTOR. Different treatments of QED corrections are chosen by 
settings of flags IMEA and IDS?. 
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Figure 10: Input organization and choice of QED treatments 
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Figure 11: The output organization 0/ HECTOR 
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4 User Options 



The following set of flags and parameters may be defined by the user in the file HECTOR. INP. 
The default values of the flags, their type, names of the corresponding variables entering the 
program with the names of host COMMON blocks one may find in appendix B, table |^. 

4.1 Selection of the scattering process 

IBEA 

IBEA=+1: anti-lepton beam (positron or 
IBEA=— 1: lepton beam (electron or fi~) 

ILEP 

ILEP=1: electron (or positron) scattering 
ILEP=2: muon scattering 

ITAR 

ITAR=1: proton target 
ITAR=2: deuteron target 

POLE 

POLE: degree of lepton beam polarization (— 1 < POLE < 1) 
EELE 

EELE: lepton beam energy in GeV 
ETAR 

ETAR: energy of target beam in GeV 
lERN 

lERN: to be used for I0PT=1 only 

IERN=0: Born cross section of a lepton scattering is defined by the user 

in the subroutine usrbrn . f 
IERN=1: default Born cross section 



4.2 Selection of measurement type and cuts 



IMEA 



=9 
=10: 



lepton measurement (leptonic variables), neutral current DIS 

jet measurement (Jaquet-Blondel variables), neutral current DIS 

jet measurement (Jaquet-Blondel variables), charged current DIS (only for IOPT=l) 

mixed variables, neutral current DIS 

double angle method, neutral current DIS (only for I0PT=1) 

lepton angle and j/j^ measurement, neutral current DIS (only for IOPT=l) 

hadron measurement, neutral current DIS 

leptonic variables, charged current DIS 

S method, neutral current DIS (only for IOPT=l) 

eS method, neutral current DIS (only for IOPT=l) 
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ZCUT 

ZCUT: 



IHCU 

IHCU=0: 
IHCU=1: 



q2CU 
Q2CU: 
HM2C 
HM2C: 

IGCU 

IGCU=0: 

IGCU=1: 



GCUT 

GCUT: 

THCl 

THCl: 

THC2 



cut on the reconstructed electron beam energy in GeV (only for I0PT=1) 
see section 12.6.1 



no cut on the hadronic final state 

rejects hadronic final states with Qj^ < Q2CU [ GeV^] and invariant hadronic mass 
(photon not included) < HM2C [GeV^]; the two desired cut values 
Q2CU and HM2C have to be set 

cut on the transferred momentum squared, Qf^, [GeV^]; see IHCU above 
cut on the hadronic invariant mass square, [GeV^]; see IHCU above 

no cut on photon variables 

rejects all final states with a photon energy > E^ut [ GeV] and a photon angle 
in the interval {9cut2, ^cuti) (angles in radians inside (0, vr) with respect to the 
electron beam); to be used only for I0PT=2, IMEA=1 

cut on the photon energy in GeV; see IGCU above 

minimal photon angle 6cut2 in rad; see IGCU above 



THC2: maximal photon angle 6cuti in rad; see IGCU above 



4.3 Selection of the main program paths and related flags 



I OPT 
IOPT= 
IOPT= 



= 1: HELIOS branch: QED LLA RC calculated in first and higher orders 
=2: TERAD branch: complete calculation of EW and QED RC in 0{a) with 
soft photon exponentiation; the model independent TERAD results may be 
supplemented by hadronic corrections from DISEP (for IDIS=2,3); 
see description of IDSP and ILOW below 
I0PT=3: combined calculation using both branches; here the TERAD results are 

supplemented by HELIOS higher order terms beginning with 0{a'^)] hadronic 
corrections from DISEP may also be added (for IDIS=2,3); the soft photon 
exponentiation of TERAD is switched off 



lORD 

IORD= 

IORD= 



IORD= 
IORD= 



Born cross section (for IOPT=l) 

QED corrections to 0{aL) (for I0PT=1) or to 0{a) with soft photon 

exponentiation (for I0PT=2) 

QED corrections to 0{{aLf) only (for I0PT=1) 

QED corrections to 0{aL) and 0{{aLY) are summed and the HELIOS 
soft photon exponentiation is included (for I0PT=1) or complete 0{a) corrections 
are summed with the (9((aL)^) and the soft photon exponentiation terms from 
HELIOS (for I0PT=3) 
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IBOS 
IBOS: 

IBOS= 
IBOS= 
IBOS= 
IBOS= 

ICOR 

ICOR: 

ICOR= 

ICOR= 

ICOR= 

ICGR= 

ICOR= 

I CMP 

ICMP: 

ICMP= 

ICMP= 

ICMP= 

lEPC 
lEPC: 
IEPC= 
IEPC= 

IWEA 

IWEA= 

IWEA= 



IVPL 




IVPL= 


=0: 


IVPL= 


= 1: 


lEXP 




IEXP= 


=0: 


IEXP= 


=1: 


IDSP 




IDSP= 


=0: 


IDSP= 


= 1: 


IDSP= 


-2: 


IDIS 




IDIS= 




IDIS= 


-2: 


IDIS= 


=3: 


IDQM 





calculation of RC to different terms in the MC cross section; may be used for all 
values of IMEA but 3,8 

only photon MC exchange contribution (|7^| term) 

only the photon - Z boson interference contribution to the MC exchange {\'^Z\ term) 
only the Z boson MC exchange contribution {\Z'^\ term) 
sum of all MC exchange contributions 

to be used for I0PT=1 only 
initial state lepton radiation 
final state lepton radiation 

Compton peak contribution (available for IMEA=1 only) 
initial plus final state radiation 

initial and final state radiation plus Compton peak contribution 

to be used for I0PT=1 only 
no Compton term 
onefold Compton integral ( |2.110| ) 
twofold Compton integral ( p.lll| ) 



to be used for I0PT=1,3 only 

e^e" conversion in 0{{aVf') is switched off 

e^e~ conversion is included 

EW RC are applied via sin^^, the effective EW mixing angle 
the full set of EW form factors is calculated; 

Caution! - using this flag may enlarge the CPU time considerably 
a = 1/137 ... is not running 

running a due to the complete fermionic vacuum polarization 
no soft photon exponentiation 

soft photon exponentiation is applied; see description of lOPT 
the DISEP branch in TERAD is not used 

together with TERAD, the DISEP branch is used (possible only for IMEA=1,4,8) 

description of flag IDIS 

only the DISEP branch is used 

DISEP leptonic corrections 

DISEP leptonic and lepton quark interference corrections 
DISEP leptonic, interference, and quarkonic corrections 



see 
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IDQM: 

IDQM= 
IDQM= 
IDQM= 

ILOW 
ILOW: 
ILOW= 
ILOW= 



to be used for IMEA=4, IDIS=2 only, see |T3 
quark masses AMQ: AMQ 
AMQ 
AMQ 



xM, M the proton mass 
const, no ISR LLA corrections from quarks 
const, no ISR or FSR LLA correction from quarks 



to be used for IMEA=1, I0PT=2 only 
no use of TERADLOW 
use of TERADLOW 



4.4 Structure function parameterizations 



IMOD 

IMOD= 

IMOD= 



=0: 
=1: 



ISTR 




ISTR: 




ISTR= 


=0 


ISTR= 


=1 


ISTR= 


=2 


ISSE 




ISSE: 




ISSE= 


= 1 


ISSE= 


=2 


ISSE= 


=3 


ISSE= 


-A 


ISSE= 


=5 


ITER 




ITER= 


=0 


ITER= 


=1 



IVAR 

IVAR<3: 



IVAR= 
IVAR= 
IVAR= 
IVAR= 

IGRP 
IGRP: 
IGRP: 



=-1 
=0: 
= 1 
=2 



use of quark distributions for the target structure functions 

the parameterization of structure functions is user supplied in the file usrstr . f 

(e.g. for a model independent analysis) 

to be used for IMOD=0 only 

use of recent parameterizations for parton distributions selected by flags ISCH 
and ISSE 

parton density functions are user supplied via file usrpdf . f 
PDFLIB is used for parton densities 



to be used for ISTR=0 
parton distributions CTEQ3 |^| 

parton distributions GRV [^, ^ (cf. the comments in sect. ( |2.3.2| )) 
parton distributions MRS (A) p5| 
parton distributions MRS (A' ) ||25| 



parton distributions MRS(G) ||25| 



quark distributions and a^^^ (if running, see flag IVPL) in leptonic bremsstrahlung 

are artifically chosen to depend on as defined from lepton momenta 

they depend on as defined from the hadronic momenta (Q^ is the integration 

variable) 

use of a damping factor for Q"^ < Q^: limQ2^o li^i Q"^) = 0; 

the value of Ql depends on the method 

the chosen quark distributions are not modified 

Volkonsky/Prokhorov damping |^ 

Stein damping p9| 

Brasse/Stein parameterization BO, 120 



to be used for ISTR=2 

a PDFLIB parameter; it defines an author group in the library (see PDFLIB 
write-up |[44|| ) 
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I SET 



T^FT- 




to hp ntiprl for T^TR — 9 


ISET: 




PDPT TR TiP{V f\vr\f^\ f^v wnipn npfinpc; f hp Qnppifir' nprpmpf priypif inn c;pf wif nin f rip (xronn 






HpfinpH hv TGRP- spp Il44in 

VJ.V!'llllv!'V-l XVJllir ^ occ II 1 1^1 y 


ISCH 






ISCH: 




to be used for ISTR=0,2 


ISCH: 




select the QCD factorization scheme; the flag has to be chosen in accordance with the 






selected parameterization in PDFACT or PDFLIB 


ISCH= 


=0 


leading order QCD 


ISCH= 


=1 


DIS factorization scheme 


ISCH= 


=2 


MS factorization scheme 


4.5 


Bins and grid 


TRIM 






TUTTVr — 
±0 ± iv — 


=0 


no integration over bins 


IBIN= 


=1 


user supplied integration over bins with subroutine usrint . f 








IBCO= 


-1 
- 1 


Ullio cLi" U-t^iiiiCU. lii Lilt; VcLildUlCo VrirVi — Jb cLliU. Viin^ — U 


TRfn— 


=2 


bins are defined in the variables VARl=x and VkK2=Q^ 


TRfn— 


o 

=o 


bins are defined in the variables VAR1=?/ and VAR2=(5^ 


ITVl 






ITV1= 


=U 


default bins in the first variable (VARl); see ([B.32|) in appendix P 


ITvl= 


=1 


loffarithmical ffrid in the first variable 


TTV1 — 


=2 


linear grid in the first variable 


TTV1 — 

X 1 V X — 


=3 


the bins in the first variable are defined by the user as a set of numbers in USRBIN; 






see VARl below 


X 1 V ± — 


-A 


user supplied bins in variable 1 via FFREAD 


ITV2 






ITV2= 


=0 


default bins in the second variable (VAR2) 


ITV2= 


=1 


logarithmical grid in the second variable 


ITV2= 


=2 


linear grid in the second variable 


ITV2= 


=3 


bins in the second variable are defined by the user as a set of numbers in USRBIN; 






see VAR2 below 


ITV2= 


-A 


user supplied bins in variable 2 via FFREAD 


NVAl 






NVAl: 




number of bins (< 100) in the first variable, to be used for ITV1=1,2 


NVA2 






NVA2: 




number of bins (< 100) in the second variable, to be used for ITV2=1,2 


VIMN 






VIMN: 




minimal value of the first variable (VIMN > 0.0), to be used for ITV1=1,2 



VIMX 
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VIMX: maximal value of the first variable (VIMX < 1.0), to be used for ITV1=1,2 
V2MN 

V2MN: minimal value of the second variable (V2MN > 0.0), to be used for ITV2=1,2 
V2MX 

V2MX: maximal value of the second variable (V2MX: y < 1.0 or < s), to be used 
for ITV2=1,2 

VARl 

VARl: user supplied bins in the first variable, a set {N < 100) of REAL numbers from 
the interval 0.0 < a; < 1.0 for IBC0=1,2, or 0.0 < y < 1.0 for IBC0=3; to be 
used for ITV1=4 

VAR2 

VAR2: user supplied bins in the second variable, a set (A^ < 100) of REAL numbers 
from the interval 0.0 < y < 1.0 for IBC0=1, or Q^^^ < < s for IBC0=2,3, 
(Qmm depends on the structure functions used, usually Q^j„=4 or 5 GeV^); 
to be used for ITV2=4 

4.6 User initialization 

lUSR 

IUSR=0: default initialization of libraries, electroweak constants, etc. 

IUSR=1: some additional initialization may be done using the external subroutine usrini.f 
where the user may initialize her/his set of wanted quantities 

4.7 Choice of an output 

I OUT 

I0UT=0,1: output is written into file HECTOR. OUT and on the screen 
I0UT=2: output is created by user subroutine usrout.f 
I0UT=3: a simple output for creation of figures 



After reading the data cards supplied by the user the program checks their compatibility. The 
following flags have highest priority: lOPT, IMOD, ITVl, and ITV2. If there are flags being in 
conflict with them they will be ignored. In other cases of incompatibility the program will 
stop and write a message into the file HECTOR . OUT. 
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4.8 User subroutines 

The user may organize her/his own interface for using the code with help of several subroutines, 
which can be rewritten by the user. Simple examples can be found inside the correspondent 
files. 

usrini.f : is called if IUSR=1. 
CALL USRINI 

In this subroutine the user may initialize parameters for later use. They can be transferred 
via COMMON /HECUSR/. 

usrdpf -f : is called if ISTR=1. 

CALL USRPDF(X,Q,UV,DV,US,DS,SS,CS,BS,TS,GLU) 

In this subroutine the user may put a set of parton density functions. 
INPUT: X,Q, with Q=v^- 
OUTPUT: UV,DV,US,DS,SS,CS,BS,TS,GLU 

Note that all density functions should be multiplied by x, i.e. UV = a;M^(a;, Q^). If the chosen 
density functions are calculated up to the next to leading order one has also to implement 
QCD radiative corrections, in dependence on the factorization scheme used. 

usrstf .f : is called if IM0D=1. 

CALL USRSTR(X,Q,FL,GL,HL,F2,G2,H2,XG3,XH3,WL,W2,XW3) 

In this subroutine the user may put a set of target structure functions. 
INPUT: X,Q, with Q=VQ^. 

OUTPUT: FL,GL,HL,F2,G2,H2,XG3,XH3,WL,W2,XW3 

Note that in this case one may not necessarily refer to the parton picture but use a model 
independent parameterization of the structure functions instead. 

usrbin.f : is called if ITV1=3 and/or ITV2=3. 
CALL USRBIN 

In this subroutine the user may define a set of points in two variables, where the cross sections 
and radiative correction are calculated. Note that the internal definitions of usrbin . f have 
higher priority than the definitions in HECTOR. INP. An example is given. 
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usrint.f : is called if IBIN=1. 
CALL USRINT 

In this subroutine the user may organize an integration over bins. A simple example is illus- 
trated by the default routine. 

usrbrn.f : is called if IBRN=0. 

CALL USRBRN(X,Y,S,SIGO) 
This subroutine can only be used in the HELIOS branch of the code (I0PT=1). 

INPUT: X,Y,S 

OUTPUT: SIGO, the Born cross section (Pa/dxdy in [nb]. 

The radiative corrections to the process defined are then calculated according to the flags 
defined in HECTOR. INP. 

usrout.f : is called if I0UT=2. 
CALL USROUT 

In this subroutine the user may organize her/his own output. A simple example, which creates 
OUTPUT VECTORS to be used e.g. in graphics programs, is given. 

usrgrf .f : is called if I0UT=3. 
CALL USRGRF 

This is another example of the user supplied OUTPUT routine creating OUTPUT VECTORS of a 
different format. 
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5 Examples of numerical output of HECTOR 



In the initialization phase, the user has to decide whether user supphed input shall be used in 
addition. If not, IUSR=0 has to be chosen (see section ^^ ). The setting of all the electroweak 
parameters is done in the file hecset.f. The unexperienced user can refer to the default 
flags. The beam parameters (particle type, charge, polarization, energy) are selected by the 
parameters given in section |4.1j . Furthermore, one has to decide in which variable set the QED 
corrections shall be calculated (flag IMEA) and whether (possible) cut parameters shall be set. 
Then, the type of corrections to be applied has to be defined, including a choice of branch(es) 
(HELIOS, TERAD, DISEP, or a combination of them) to be used. This is done by the flags 
I OPT, lORD etc. being described in section [4.3| . The structure function parameterizations are 
chosen by the flags of section [4.4| . For IBIN=0, the output is calculated for the pre-defined 
grid of kinematical points and no integration over bins is performed. The flags and parameters 
of section E75| have to be set. 



Two example files for HECTOR. INP, HECTOR . INP . SRT and HECTOR. I NP.STD are available. 
The sample outputs (appendices C, D) correspond to these input files (renamed as HECTOR. INP). 
In addition, we show in tables 3-6 cross section values for varying input parameters and in 
figures p!^-[14| radiative corrections in percent for three selected types of variables. These 
examples are foreseen to act as benchmarks for typical variations in the use of HECTOR. 



ISSE 


ISCH 


cxb [nb] 


6^'^ 

"mi 


[%] 




1 





0.2246E+06 


13.01 


2 





0.2441E+06 


12.94 


1 


1 


0.2225E+06 


12.65 


3 


1 


0.2679E+06 


11.76 


1 


2 


0.2345E+06 


12 


.82 




2 


2 


0.2702E+06 


12.13 


3 


2 


0.2695E+06 


11 


.84 


4 


2 


0.2410E+06 


13.56 


5 


2 


0.2415E+06 


13.55 



Table 3: Born cross sections for unpolarized e^p scattering with NLO QCD corrections for 
different sets of parton parameterizations. Parameters: x = 10~^,y = 0.1, = 26.7 GeV, 
Ep = 820.0 GeV, and IVPL = 1 (i.e. running a). The low damping corresponds to IVAR 
= 0. 



IVAR 


(Tb [nb] 


6^'^ 


[%] 


-1 


0.4026E+08 


9.452 





0.3816E+08 


9.222 


1 


0.3616E+08 


9.197 


2 


0.7445E+07 


32.53 



Table 4: Born cross sections for e~^p scattering with different parton distribution parameter- 
izations at low Q^. Parameters: x = 10'^, y = 0.1, E^ = 26.7 GeV, Ep = 820.0 GeV, ISSE 
= 1, ISCH = 0, IVPL = 1. 
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lOPT 


IDSP 


IWEA 


as [nb] 




[%] 


3 
3 







1 


0.1260E-02 
0.1278E-02 


-3.980 
-5.130 


2 
2 


2 
2 



1 


0.1260E-02 
0.1278E-02 


-5.129 
-6.268 



Table 5: Born cross sections for e^p scattering for different values of flag IWEA. Parameters: 
X = 0.5, y = 0.5, = 26.7 GeV, Ep = 820.0 GeV, ISSE = 1, ISCH = 0, IVAR = 0, IVPL = 1. 



IMEA 


ICOR 


IDSP 


I OPT 


cxb [nb] 


5(1) 


[%] 




[%] 




[%] 


Comments 


1 


1 





1 


0.2246E+06 


13.01 


-.4092 


12.60 




2 


1 





1 


0.2246E+06 


3.928 


.3942 


4.322 




3 


1 





1 


0.1066E+01 




m 


.4420 


.3583 




4 


1 





1 


0.2246E+06 


10.22 


.6968 


10.92 




5 


1 





1 


0.2246E+06 


5.512 


.2490 


5.761 


E^^t = 35 GeV 


6 


1 





1 


0.2246E+06 


6.182 


.2334 


6.416 


E,,.t = 35 GeV 


7 


1 





1 


0.2246E+06 


4.469 


.4151 




^84 




8 


1 





1 


0.1066E+01 


3.172 


-.1664 


3.006 




9 


1 





1 


0.2246E+06 


4.150 


.5021 


4.652 




10 


1 





1 


0.2246E+06 


9.328 


2.492 


11 


.82 




1 


2 





1 


0.2246E+06 


4.643 












4 


2 





1 


0.2246E+06 


-2.687 












9 


2 





1 


0.2246E+06 




m 












10 


2 





1 


0.2246E+06 


-3.040 












1 


3 





1 


0.2246E+06 


.02006 










ICMP=1 










0.2246E+06 


.02753 










ICMP=2 


1 


1 





2,3 


0.2246E+06 


17.82 






17.41 




2 


1 





2,3 


0.2246E+06 


3.748 






4.133 




4 


1 





2,3 


0.2246E+06 


7.195 






7.783 




7 


1 





2,3 


0.2246E+06 


4.337 






4.736 




1 


1 


1 


2,3 


0.2246E+06 


17.82 






17.41 




1 


1 


2 


2,3 


0.2246E+06 


17.95 






17.41 




4 


1 


1 


2,3 


0.2246E+06 


7.189 






7.783 




4 


1 


2 


2,3 


0.2246E+06 


7.189 






7.783 




8 


1 


1 


2,3 


0.1066E+01 


8.100 






7.933 


IDIS=2 










0.1066E+01 


8.070 






7.903 


IDIS=1 



Table 6: Cross sections and radiative corrections for e~^p scattering as functions of flags IMEA, 
ICOR, IDSP, lOPT. Parameters: x = 10~^?/ = 0.1, Ee = 26.7 GeV, Ep = 820.0 GeV, and IVPL 
= 1 (i.e. running a). The parton density parameterization CTEQ3L0 is chosen together with 
the low damping factor corresponding to IVAR = 0. Longitudinal structure functions are 
disregarded. For the 0(a^) corrections e+e~ beam conversion terms were added. For lOPT = 2 
(3), the values of 5^^^ are given in columns 6(8). Note the comments in 1^ for the branches 
selected by IMEA = 5,6. 
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Figure 12: Radiative corrections in leptonic variables for the default settings in percent. Dotted 
lines: 0{a), dashed lines: 0{a^), solid lines: in addition soft photon exponentiation. 
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Figure 13: Radiative corrections in mixed variables for the default settings in percent. Dotted 
lines: 0{a), dashed lines: 0{p?\ solid lines: in addition soft photon exponentiation. 
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Figure 14: Radiative corrections in double angle variables for the default settings in percent. 
Dotted lines: 0{a), dashed lines: 0{a^), solid lines: in addition soft photon exponentiation. 
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6 Output 



The code writes a file HECTOR. OUT by default (I0UT=0). The user can organize her/his own 
output in the user file usrout .f , which is called for I0UT=2. Option I0UT=3 pipelines output 
to the user file usrgrf .f. Examples for these subroutines and some instructions on their use 
are given in the files usrout . f and usrgrf . f . Examples of the program output may be found 
in appendices C, D. 

The program writes the following messages into the file HECTOR . OUT: 

• program Logo labeled with the version number; 

• the data cards which have been read in the FFREAD format; 

• in case of an incompatible fiag definition the program writes the values of 
the corresponding fiags into the file HECTOR . OUT and terminates thereafter; 

• version number of the library PDFLIB (if activated by the user); 

• an information about the chosen program branch; 

• type of process and variables selected; 

• target type; 

• information about the projectile particle: its charge, type and polarization; 

• type the structure functions used; 

• choice of the low Q'^ modification used; 

• input parameters and calculated electroweak parameters; 

• energy characteristics of the colliding particles; 

• ZCUT cut parameter (for I0PT=1); 

• choice of bin variables VARl and VAR2; 

• type of grid in each variable; 

• number of bins in each variable; 

• label of the chosen set of parton densities (for ISTR=0). 

The results of the calculation are presented in six columns: values of the first variable VARl, 
values of the second variable VAR2, the s values in GeV^, Born cross section values in nbarn, 
corrected cross section values in nbarn, correction values in %. 
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A Definition of scaling variables in presence of hard 
photon emission 

In this appendix, we collect the definitions of all the sets of kinematical variables, which may 
be used in HECTOR. Some motivations for the specific choices may be found in [H-fH. 

Leptonic variables: 

The scaling variables are defined from leptonic variables: 

yi = --I , (A.l) 

Q2 = (k, - = -2k,k2, (A.2) 

XI = (A.3) 
Syi 

Hadronic variables: 

The scaling variables are defined from hadronic variables using the 4-vector of the current jet: 

Vh = ; , (A.4) 

Ql = (P2-Pl)^ (A.5) 
XH = (A.6) 



Jaquet-Blondel variables 

Again, the scaling variables are defined from hadronic variables. Using the transverse hadronic 
momentum in the laboratory system, p^, instead of p2 and the energy fiow in z direction in 
the definition of the scaling variables^ : 



Q% = (A.8) 
^.B = (A.9) 

with 

E = ^(E;,-p;,,,), (A.IO) 

h 



^Note that in the complete 0{a) calculations for Jaquet-Blondel and mixed variables in refs. ||T^] the 
approximate relation tJj^ « yh was used which holds iff all masses in the final state can be neglected with 
respect to the particle momenta. 
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and Eh,pfi^x,ph^y,ph^2 being the four-momentum components of the hadronic final state. 

Mixed variables /f^/; 

The scahng variables are defined from both leptonic and Jaquet-Blondel variables: 



(A.ll) 
(A.12) 

(A. 13) 



For mixed variables, the allowed values of x are not necessarily restricted to be smaller than 



Qm 


= Qh 


Vm 


= y.m^ 




Qf 


2'm 





one; see for details appendix B.3.1 of |]TT 



E method 



The scaling variables are defined from both leptonic and hadronic variables as being measured 
in the laboratory system: 



J: + E'l {1- COS Oi)' 



Ql 



Sy^ 



(A.14) 
(A.15) 
(A.16) 



The integration boundary in the case of final state radiation is the solution of the equation 



eI^^ + E2E1 



-1/3 



with 



El 
El 

E\ 
E2 



x{El + E\] 



--y\^--y-^xy J + -xy [I- y), 



2 

1 

6 
2 

-xy [l-y+ -xy^ 



'9(1 + +2?/2 
1 



27 - 18y - '^xy^ + 2xy Q - y + y^ 



(A. 17) 
(A.18) 

(A.19) 
(A.20) 

(A.21) 
(A.22) 
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eH method j^/; 

The eS method mixes variables from the set of leptonic variables and of the S method: 

C = Ql (A-23) 
= a;^, (A.24) 

= ^ = 1^- ^^-'^^ 
The integration boundary Zq'-^ is the same as in the S method. 

Double angle method /j^/; 

The scaling variables are defined from both leptonic and hadronic scattering angles as being 
measured in the laboratory system. 

The definitions of scaling variables in the double angle method are: 

4i?fcos2| 

sm h sm — cos — tan — 

2 2 2 2 

01 



sm 
2 



01 61 9h ' 
sm h cos — tan — 

2 2 2 



(A.27) 



= (A.28) 

6^2/ method 

A quite similar opportunity is to express the cross sections in terms of 6e in the laboratory 
system and The kinematical variables are: 



= AEtil-yJ- -i, (A.29) 

' 1 — cos Oi 

yoj = VjB^ (A.30) 
xey = (A.31) 

^Udy 
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Table of user options 





FFRFAD flap- 




Tiarampf PT 


pom m on 


f] pfa n 1 f 




1 
i 




T) 

tx 




niLCr AK 


oU.U 


A 1 
4. i 


9 


T7TAD 

r,i An 


r> 
IX 


f DD 


nn,L.r AK 


oZU.U 


A 1 
4.1 


o 
o 


nnTTT 


R 
IX 






u.u 


/I 
4.Z 


A 
4 




r> 
IX 


A MIT OPT" 


nrLUnUL. 


n n 
U.U 


A 
4.Z 





iDL.U 




TDPn 
ibUU 




i 


/I K 
4.D 





TTDT? A 




TOT7 A M 

IdEjAM 


nr,UUr i 


1 

— i 


/I 1 
4. i 


7 






TDTAT 

IhSliN 


IJT7PDTAT 


n 
U 


A c; 
4.0 


o 
o 


TDnc 
IdUo 




TTT7DM 




4 


/I Q 
4.0 


Q 


TT3DM 




TTTT3BM 




u 


/I Q 
4.0 


1 n 

iU 






IL-rlr 1 




i 


A Q 
4.0 


i i 


TPnD 
lUUri 











/I Q 
4.C) 




Tr\TC! 

IJJlD 






UHTPTnC! 


i 


/I Q 

4.0 












i 


/I Q 

4.0 


1 /I 


lUbr 




lUbr 


niiL-L.UK 


u 


/I Q 

4.0 


iO 










n 
U 


A Q 

4.0 


iO 


TT7YD 




TT7YD 




n 
U 


/I Q 
4.C) 


1 7 


TPr'TT 




IbUU 1 




n 
U 


/I o 
4.Z 


io 


i(jKr 




JMbnUUr 


risLLitrUr 


6 


/I A 
4.4 


1 O 


InUU 




InOU 1 


wunTur* 


U 


/I O 
4.Z 


on 


T T ITD 




iLbr i U 


nhOUr 1 


1 


1 1 
i. i 


oi 
zi 


TT ni.T 
ILUW 




TT n\j 


niiULUW 


U 


A Q 

4.0 


zz 


TMT? A 

InrLR 




TMT7 A O 


nhUUr i 


i 


/I O 
4.Z 


Zo 


IrlUU 




rlUUr,L 


Uf PT7T P 


n 
U 


/I A 
4.4 


0/1 
Z4 


T DDT 

iUr i 




lUr 1 


nilCUr 1 


i 


/I Q 
4.C) 


O K 
ZO 


iUnl) 




T nDr\ 


uu'ppnD 
rlEOUUK 


i 


A Q 

4.0 


Ofi 
ZD 


lUU i 




THTTT 


UTTPnTTD 

llllL.UUr 


u 


/I 7 
4. / 


27 


ISCH 




ISCHM 


HECQCD 


2 


4.4 


28 


I SET 




NSET 


HECPDF 


41 


4.4 


29 


ISSE 




ISSET 


HECSTF 


1 


4.4 


30 


ISTR 




ISTR 


HECSTF 





4.4 


31 


ITAR 




ITARG 


HECOPT 


1 


4.1 


32 


ITER 




ITERAD 


HECCOR 


1 


4.4 


33 


ITVl 




ITVl 


HECGRD 





4.5 



Table 7: C/ser defined flags in the file HECTOR. INF 
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rrKLAU nag 


ty])0 


])araiiiotor 


COllllllOll 


,1 ,,r , , ,1-1- 

(Ic'iaiiit 


roioroiic'cs 


34 


ITV2 




ITV2 


HECGRD 


2 


4.5 


35 


lUSR 




lUSR 


HECUSR 





4.6 


36 


IVAR 




IVAR 


HECFLG 





4.4 


37 


IVPL 




IVPOL 


HECVPL 


1 


4.3 


38 


IWEA 




IWEAK 


HECFLG 





4.3 


39 


NVAl 




NVl 


HECBIN 


6 


4.5 


40 


NVA2 




NV2 


HECBIN 


20 


4.5 


41 


POLE 


R 


POL 


HECPAR 


0.0 


4.1 


42 


THCl 


R 


THCUTl 


HECGMC 


0.0 


4.2 


43 


THC2 


R 


THCUT2 


HECGMC 


0.0 


4.2 


44 


VARl 


R 


Vl(lOO) 


HECBIN 




4.5 


45 


VAR2 


R 


V2(100) 


HECBIN 




4.5 


46 


VIMN 


R 


VIMI 


HECGRD 


l.OE-04 


4.5 


A 7 

4 1 


VIMa 


ix 


V IMA 




i.U 


4.0 


48 


V2MN 


R 


V2MI 


HECGRD 


l.OE-02 


4.5 


49 


V2MX 


R 


V2MA 


HECGRD 


1.0 


4.5 


50 


Q2CU 


R 


TCUT 


HECHDC 


0.0 


1.2 


51 


ZCUT 


R 


EECU 


HECPAR 


35.0 


4.2 



Table 7 (continued): The flags to be set by the user in the file HECTOR. INP. 



The default bins in x are: 

X = 10-^ 10-^ 10"^ 10-\ 0.5, 0.9 . (B.32) 

The default values of y are defined in table 7. They are located at 20 equidistant points in 
the interval from 0.01 to 1. 
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C Example file HECTOR . OUT 

The following sample file HECTOR. OUT corresponds to the first row of Table 3. 

HECTOR vl.OO 

1 USER'S DIRECTIVES TO RUN THIS JOB 





r> ATA 


r* ARn 




T TCIT 






riATA 


TARn 




Lj . . . . 


QFTTT 
OEil 1 ±. 




UR L ii 


r" ARn 


f niVTTPMT 








n AT A 


TARn 
v^iinjj 


f nWTFMT 


n 

\j . , . , 


r J_jiiUO 




DATA 


TARn 
v-/iinjj 




THPT 






r> ATA 
Uiilii 


r* ARn 




TMP A 
IrlEiii 


1 




r> ATA 
Uiiiii 


r* ARn 




TRF A 
±DEiii 


1 
1 




nATA 
JJii i ii 


fARn 
L'iiriJJ 


f niVTTPMT 


TT PP 
J.J_jEjr 


1 




n AT A 

JJii i ii 


TARn 
v^iinjj 


f nWTFMT 


TTAR 

± 1 iin 






nATA 
UR. i ii 


r* ARn 


f nWTPWT 








r> ATA 


r" ARn 




THRn 


1 
1 




riATA 

UR 1 R 


TARn 






A 
*± 




nATA 

JJii i ii 


fARn 
v^iinjj 


f nWTFMT 


Tf HR 


1 




nATA 

JJii i ii 


fARn 




Tf MP 

J. L*! Jr 


1 




nATA 

Uii 1 ii 


fARn 


f DMTFWT 


TPPf 






nATA 


r* ARn 


f nMTT7MT 
OUIM 1 EjIM 1 


TRRM 


■1 
1 




nATA 

Uii i ii 


fARn 


OUIM 1 EjIM 1 








nATA 
JJii i ii 


fARn 
UiiriJJ 


f nWTFMT 


IJJor 






nATA 

JJii i ii 


fARn 
v^iinjj 


V^UJM i EjJM i 


Tnic; 

J-U J-O 






nATA 

UR 1 ii 


fARn 
uiinjj 


f niUTPWT 

LfUJ)! 1 EjIM 1 


TnnM 

J.Ul^J J 






DATA 


CARD 


CONTENT 


ILOW 







DATA 


CARD 


CONTENT 








DATA 


CARD 


CONTENT 


IMOD 







DATA 


CARD 


CONTENT 


ISTR 







DATA 


CARD 


CONTENT 


ISSE 


1 




DATA 


CARD 


CONTENT 


ITER 


1 




DATA 


CARD 


CONTENT 


IVAR 







DATA 


CARD 


CONTENT 


IGRP 


3 




DATA 


CARD 


CONTENT 


I SET 


17 




DATA 


CARD 


CONTENT 


ISCH 







DATA 


CARD 


CONTENT 








DATA 


CARD 


CONTENT 


IWEA 







DATA 


CARD 


CONTENT 


IVPL 


1 




DATA 


CARD 


CONTENT 


lEXP 







DATA 


CARD 


CONTENT 








DATA 


CARD 


CONTENT 


IHCU 
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DATA 


CARD 


CONTENT 


Q2CU 


0. 





***** 


DATA 


CARD 


CONTENT 


HM2C 


0. 





***** 


DATA 


CARD 


CONTENT 








***** 


DATA 


CARD 


CONTENT 


IGCU 







***** 


DATA 


CARD 


CONTENT 


GCUT 


0. 





***** 


DATA 


CARD 


CONTENT 


THCl 


0. 





***** 


DATA 


CARD 


CONTENT 


THC2 


0. 





***** 


DATA 


CARD 


CONTENT 








***** 


DATA 


CARD 


CONTENT 


C 


PARAMETERS 


***** 


DATA 


CARD 


CONTENT 


POLE 


0. 





***** 


DATA 


CARD 


CONTENT 


EELE 


2e 


1.70 


***** 


DATA 


CARD 


CONTENT 


ETAR 


820.0 


***** 


DATA 


CARD 


CONTENT 


ZCUT 


C 


1.00 


***** 


DATA 


CARD 


CONTENT 








***** 


DATA 


CARD 


CONTENT 


C 


LATTICE 


***** 


DATA 


CARD 


CONTENT 


IBIN 







***** 


DATA 


CARD 


CONTENT 


IBCO 


1 




***** 


DATA 


CARD 


CONTENT 


ITVl 


4 




***** 


DATA 


CARD 


CONTENT 


ITV2 


4 




***** 


DATA 


CARD 


CONTENT 


NVAl 


1 




***** 


DATA 


CARD 


CONTENT 


NVA2 


1 




***** 


DATA 


CARD 


CONTENT 


VIMN 


1. 


E-4 


***** 


DATA 


CARD 


CONTENT 


VIMX 


1. 





***** 


DATA 


CARD 


CONTENT 


V2MN 


0. 


01 


***** 


DATA 


CARD 


CONTENT 


V2MX 


0. 


5 


***** 


DATA 


CARD 


CONTENT 


VARl 


0. 


001 


***** 


DATA 


CARD 


CONTENT 


VAR2 


0. 


1 


***** 


DATA 


CARD 


CONTENT 








***** 


DATA 


CARD 


CONTENT 


I OUT 









DATA 


CARD 


CONTENT 










DATA 


CARD 


CONTENT 


lUSR 









DATA 


CARD 


CONTENT 










DATA 


CARD 


CONTENT 


STOP 







BINS: 

VARK 1)= .1000E-02 
VAR2( 1)= .lOOOE+OO 

***** LLA & NTLA BRANCH 

***** 0(ALPHA*L) CORRECTIONS IN LLA 

***** LEPTON MEASUREMENT, NEUTRAL CURRENT 

***** CHARGED ANTI-LEPTON SCATTERING 

***** PROTON TARGET 

***** LEPTON BEAM POLARIZATION= .0 

***** L-BEAM: E 

***** PARTON DESCRIPTION FOR STRUCTURE FUNCTIONS IS USED 
***** CTEQ3 PDF PARAMETERIZATION, LO 

***** LOW Q2 MODIFICATION: Q2 < Q2MIN: VOLKONSKY/PROKHOROV DAMPING 
***** q2MIN= 4.000 GeV**2 
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***** INPUT ELECTOROWEAK PARAMETERS USDED: 
***** z-BOSON MASS AMZ (GEV) = 91.175 

***** W-BOSON MASS AMW (GEV) = 80.399 
***** HIGSS BOSON MASS AMH (GEV) = 300.000 
***** TOP QUARK MASS AMT (GEV) = 180.000 
***** CALCULATED EW PARAMETERS: 
***** SIN(THETA_W) = .232000 
***** DELTA.R = .036948 

***** BEAM KINEMATICS ***** 

***** LEPTON BEAM ENERGY (GEV) = 26.70000076293945 

***** PROTON BEAM ENERGY (GEV) = 820.0 

***** BINNING IN X & Y 

***** USER SUPPLIED BINS IN VARIABLE 1 

***** USER SUPPLIED BINS IN VARIABLE 2 

***** NUMBER OF BINS IN VI: NVA1= 1 

***** NUMBER OF BINS IN V2: NVA2= 1 

CTEQ3(L0 ) 

VARl VAR2 S BORN SIG.CORR DELTA 

.lOOE-02 .lOOE+00 .876E+05 .225E+06 .254E+06 . 130E+02 
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D Test run output 

HECTOR vl.OO 



USER'S DIRECTIVES TO RUN THIS JOB 





DATA 


CARD 


CONTENT 


LIST 




***** 


DATA 


CARD 


CONTENT 


C. . . . 


SETTi: 


***** 


DATA 


CARD 


CONTENT 






***** 


DATA 


CARD 


CONTENT 


C. . . . 


FLAGS 


***** 


DATA 


CARD 


CONTENT 


I OPT 


3 


***** 


DATA 


CARD 


CONTENT 


IMEA 


1 


***** 


DATA 


CARD 


CONTENT 


IDEA 


1 


***** 


DATA 


CARD 


CONTENT 


ILEP 


1 


***** 


DATA 


CARD 


CONTENT 


I TAR 


1 


***** 


DATA 


CARD 


CONTENT 






***** 


DATA 


CARD 


CONTENT 


lORD 


3 


***** 


DATA 


CARD 


CONTENT 


IBOS 


4 


***** 


DATA 


CARD 


CONTENT 


ICOR 


5 


***** 


DATA 


CARD 


CONTENT 


I CMP 


1 


***** 


DATA 


CARD 


CONTENT 


lEPC 





***** 


DATA 


CARD 


CONTENT 


IBRN 


1 


***** 


DATA 


CARD 


CONTENT 






***** 


DATA 


CARD 


CONTENT 


IDSP 





***** 


DATA 


CARD 


CONTENT 

^/^JXV X X-lX« X 


IDIS 


1 


***** 


DATA 


CARD 


CONTENT 


IDQM 


1 


***** 


DATA 


CARD 


CONTENT 


ILOW 





***** 


DATA 


CARD 


CONTENT 






***** 


DATA 


CARD 


CONTENT 


IMOD 





***** 


DATA 


CARD 


CONTENT 

^/^JXV X X-lX« X 


ISTR 





***** 


DATA 


CARD 


CONTENT 


ISSE 


1 


***** 


DATA 


CARD 


CONTENT 


ITER 


1 


***** 


DATA 


CARD 


CONTENT 


IVAR 





***** 


DATA 


CARD 


CONTENT 


IGRP 


3 




n AT A 


r* ARD 






1 7 




DATA 


CARD 


CONTENT 


ISCH 







DATA 


CARD 


CONTENT 








DATA 


CARD 


CONTENT 


IWEA 







DATA 


CARD 


CONTENT 


IVPL 


1 




DATA 


CARD 


CONTENT 


lEXP 


1 




DATA 


CARD 


CONTENT 








DATA 


CARD 


CONTENT 


IHCU 







DATA 


CARD 


CONTENT 


02CU 


0.0 




DATA 


CARD 


CONTENT 


HM2C 


0.0 




DATA 


CARD 


CONTENT 








DATA 


CARD 


CONTENT 


IGCU 
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***** DATA CARD CONTENT GCUT 0.0 
***** DATA CARD CONTENT THCl 0.0 
***** DATA CARD CONTENT THC2 0.0 
***** DATA CARD CONTENT 

***** DATA CARD CONTENT C PARAMETERS 

***** DATA CARD CONTENT POLE 0.0 
***** DATA CARD CONTENT EELE 26.70 
***** DATA CARD CONTENT ETAR 820.0 
***** DATA CARD CONTENT ZCUT 0.00 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 

***** DATA CARD CONTENT VARl 0.0001 0.01 0.1 0.5 
***** DATA CARD CONTENT VAR2 0.01 0.1 0.5 0.9 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT 
***** DATA CARD CONTENT STOP 
BINS: 

VARK 1)= O.lOOOE-03 
VARK 2)= O.lOOOE-01 
VARK 3)= O.lOOOE+00 
VARK 4)= 0.5000E+00 
VARK 5)= 0.9000E+00 
VAR2( 1)= O.lOOOE-01 
VAR2( 2)= O.lOOOE+00 
VAR2( 3)= 0.5000E+00 
VAR2( 4)= 0.9000E+00 
VAR2( 5)= 0.9900E+00 

***** (ALPHA) COMPLETE CORRECTIONS SUPPLEMENTED 
BY LLA TERMS IN HIGHER ORDERS 

***** LEPTON MEASUREMENT, NEUTRAL CURRENT 

***** CHARGED ANTI-LEPTON SCATTERING 

***** PROTON TARGET 

***** LEPTON BEAM POLARIZATION= . OOOOOOOOOOOOOOOOE+00 
***** L-BEAM: E 



c 


LATTICE 


IBIN 







IBCO 


1 




ITVl 


4 




ITV2 


4 




NVAl 


1 




NVA2 


1 




VIMN 


1. 


,E-4 


VIMX 


1. 


,0 


V2MN 


0. 


,01 


V2MX 


0, 


,5 


VARl 


0, 


,0001 


VAR2 


0, 


,01 


I OUT 


1 




lUSR 








0.90 
0.99 
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PARTON DESCRIPTION FOR STRUCTURE FUNCTIONS IS USED 
CTEQ3 PDF PARAMETERIZATION, LO 

LOW Q2 MODIFICATION: Q2 < Q2MIN: VOLKONSKY/PROKHOROV DAMPING 
Q2MIN= 4.000 GeV**2 
INPUT ELECTOROWEAK PARAMETERS USDED: 
Z-BOSON MASS AMZ (GEV) = 91.175 
W-BOSON MASS AMW (GEV) = 80.399 
HIGSS BOSON MASS AMH (GEV) = 300.000 
TOP QUARK MASS AMT (GEV) = 180.000 
CALCULATED EW PARAMETERS: 
SIN(THETA_W) = 0.232000 
DELTA.R = . 036948 

BEAM KINEMATICS ***** 

LEPTON BEAM ENERGY (GEV) = 26.70000076293945 
PROTON BEAM ENERGY (GEV) = 820.0000000000000 
BINNING IN X & Y 
USER SUPPLIED BINS IN VARIABLE 1 
USER SUPPLIED BINS IN VARIABLE 2 
NUMBER OF BINS IN VI: NVA1= 5 
NUMBER OF BINS IN V2: NVA2= 5 



CTEQ3(L0 ) 
VARl 
O.lOOE-03 
O.lOOE-03 
O.lOOE-03 
. lOOE-03 
O.lOOE-03 
O.lOOE-01 
O.lOOE-01 
O.lOOE-01 
O.lOOE-01 
O.lOOE-01 
O.lOOE+00 
O.lOOE+00 
O.lOOE+00 
. lOOE+00 
O.lOOE+00 
. 500E+00 
0.500E+00 
. 500E+00 
. 500E+00 
. 500E+00 
0.900E+00 
. 900E+00 
. 900E+00 
. 900E+00 
. 900E+00 



VAR2 
O.lOOE-01 
O.lOOE+00 
. 500E+00 
0.900E+00 
. 990E+00 
O.lOOE-01 
O.lOOE+00 
0.500E+00 
0.900E+00 
0.990E+00 
O.lOOE-01 
O.lOOE+00 
0.500E+00 
0.900E+00 
0.990E+00 
O.lOOE-01 
O.lOOE+00 
. 500E+00 
0.900E+00 
0.990E+00 
O.lOOE-01 
O.lOOE+00 
0.500E+00 
. 900E+00 
. 990E+00 



S 

0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
. 876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 
0.876E+05 





BORN 


0, 


. 113E+10 


0, 


.382E+08 


0, 


, 113E+07 


0. 


, 324E+06 


0. 


,271E+06 


0, 


. 137E+06 


0. 


. 175E+04 


0, 


, 573E+02 


0. 


, 150E+02 


0. 


, 123E+02 


0, 


. 120E+04 


0, 


. llOE+02 


0, 


, 283E+00 


0. 


,611E-01 


0. 


,494E-01 


0, 


.956E+01 


0, 


.717E-01 


0, 


, 126E-02 


0. 


, 149E-03 


0, 


, 112E-03 


0, 


.661E-02 


0, 


.387E-04 


0. 


.546E-06 


0. 


, 540E-07 


0. 


,400E-07 



SIG.CORR 
0. 109E+10 
. 429E+08 
. 153E+07 
. 684E+06 
0.176E+07 
. 134E+06 
. 200E+04 
0.785E+02 
. 343E+02 
. 229E+03 
0.941E+03 
0. 109E+02 
. 348E+00 
0.115E+00 
0.594E+00 
0.514E+01 
0.524E-01 
0.121E-02 
0.200E-03 
0.471E-03 
. 204E-02 
. 183E-04 
0.362E-06 
0.493E-07 
0.523E-07 



DELTA 
-0.353E+01 
0. 124E+02 
. 348E+02 
O.lllE+03 
0.549E+03 
-0.233E+01 
. 139E+02 
. 369E+02 
0.129E+03 
0.175E+04 
-0.216E+02 
-0.971E+00 
. 232E+02 
. 878E+02 
O.llOE+04 
-0.463E+02 
-0.269E+02 
-0.398E+01 
. 344E+02 
0.319E+03 
-0.691E+02 
-0 . 528E+02 
-0 . 337E+02 
-0.857E+01 
0.307E+02 
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